ACT-03-03 
FTUV-030826 
KCL-PH/TH 030826 
MIFP-03-19 
UA/NPPS 26-08-03 

WMAPing the Universe: Supersymmetry, Dark Matter, Dark Energy, 
Proton Decay and Collider Physics 

A.B. Lahanas", N.E. Mavromatos^'^ and D.V. Nanopoulos"''^'^ 

" University of Athens, Physics Department, Nuclear and Particle Physics section, 
Panepistimioupolis, Zografou Campus, Athens GR-15771, Greece. 
^ Kings College London, University of London, Department of Physics, Strand WC2R 

2LS, London, U.K. 

Departamento de Fisica Teorica, Universidad de Valencia, E-46100, Burjassot, 

Valencia, Spain. 

Ceorge P. and Cynthia W. Mitchell Institute of fundamental Physics, Texas A & M 

University, College Station, TX 77843-4242, USA. 
^ Astroparticle Physics Group, Houston Advanced Research Center (HARC), The 
Mitchell Campus, The Woodlands, TX 77381, USA 
^ Chair of Theoretical Physics, Academy of Athens, Division of Natural Sciences, 28 
Panepistimiou Avenue, Athens GR-10679, Greece. 

Abstract 

In this review we critically discuss constraints on minimal supersymmetric models of particle 
physics as implied by the recent astrophysical observations of WMAP satellite experiment. Al- 
though the prospects of detecting supersymmetry increase dramatically, at least within the context 
of the minimal models, and 90% of the available parameter space can safely be reached by the 
sensitivity of future colliders, such as Tevatron, LHC and linear colliders, nevertheless we pay 
particular emphasis on discussing regions of the appropriate phase diagrams, which -if realized in 
nature- would imply that detection of supersymmetry, at least in the context of minimal models, 
could be out of colliders reach. We also discuss the importance of a precise determination of the 
radiative corrections to the muon anomalous magnetic moment, (7^ — 2, both theoretically and 
experimentally, which could lead to elimination of such "out of reach" regions in case of a con- 
firmed discrepancy oi — 2 from the standard model value. Finally, we briefly commend upon 
recent evidence, supported by observations, on a dark energy component of the Universe, of as yet 
unknown origin, covering 73% of its energy content. To be specific, we discuss how supergravity 
quintessence (relaxation) models can be made consistent with recent observations, which may lead 
to phenomenologically correct constrained supersymmetric models, accounting properly for this 
dark energy component. We also outline their unresolved problems. 



1 INTRODUCTION 



Particle Physics and Astrophysics have been separate for a number of years. Until a few 
years ago the accuracy with which astrophysical measurements were made was much lower 
than the corresponding one in particle physics experiments, thereby preventing a fruitful 
interaction between the two communities. 

However, in the past decade we have witnessed spectacular progress in precision mea- 
surements in astrophysics as a result of significant improvements in terrestrial and ex- 
traterrestrial instrumentation. The (second phase of the) Hubble telescope opened up 
novel paths in our quest for understanding the Universe, by allowing observations on 
distant corners of the observable Universe that were not accessible before. 

From the point of view of interest to particle physics, the most spectacular claims 
from astrophysics came five years ago from the study of distant supernovae (redshifts 
z ~ 1) by two independent groups [l]. These observations pointed towards a current 
era acceleration of our Universe, something that could be explained either by a non-zero 
cosmological constant in a Friedman-Robertson- Walker-Einstein Universe, or in general 
by a non-zero dark energy component, which could even be relaxing to zero (the data are 
consistent with this possibility). This claim, if true, could revolutionize our understanding 
of the basic physics governing fundamental interactions in Nature. Indeed, only a few 
years ago, particle theorists were trying to identify (alas in vain!) an exact symmetry of 
nature that could set the cosmological constant (or more generally the vacuum energy) to 
zero. Now, astrophysical observations point to the contrary. 

The skeptics may question the accuracy of the supernovae observations, given that 
neither the nuclear physics associated with their evolution, nor the physics involved in 
the intergalactic and interstellar matter are well understood to date so as to exclude the 
possibility that the observed effects on the dimering of the distant supernovae (z ~ 1) 
as compared to the nearby ones {z <C 1) are due to conventional physics and are not 
related to the geometry of the Universe. However, there is additional evidence from quite 
different in origin astrophysical observations, those associated with the measurement of 
the cosmic microwave background radiation (CMB), which point towards the fact that 
73 % of the Universe vacuum energy consists of a dark (unknown) energy substance, in 
agreement with the (preliminary) supernovae observations. Moreover, recently j2j two 
more distant supernovae have been discovered {z > 1), exhibiting similar features as the 
previous measurements, thereby supporting the geometric interpretation on the accelera- 
tion of the Universe today, and arguing against the nuclear physics or intergalactic dust 
effects. 

Above all, however, there are the very recent data from a new probe of Cosmic 
Microwave Background Radiation Anisotropy (Wilkinson Microwave Anisotropy Probe 
(WMAP)) (2. In its first year of running WMAP measured CMB anisotropics to an 
unprecedented accuracy of billionth of a Kelvin degree, thereby correcting previous mea- 
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surements by the Cosmic Background Explorer (COBE) satellite jj] by several orders of 
magnitude. This new satellite experiment, therefore, opened up a new era for astroparticle 
physics, given that such accuracies allow for a determination (using best fit models of the 
Universe) of cosmological parameters j5|, and in particular cosmological densities, which, 
as we shall discuss in this review, is quite relevant for constraining models of particle 
physics to a significant degree. 




Figure 1: The energy content of our Universe as obtained by fitting data of WMAP satel- 
lite. The chart is in perfect agreement with earlier claims made by direct measurements 
of a current era acceleration of the Universe from distant supernovae type la (courtesy of 
\http://map.gsfc.nasa.gov/\} . 

The WMAP satellite experiment determined the most important cosmological param- 
eters that could be of relevance to particle physicists, namely [3j: the Hubble constant, 
and thus the age of the Universe, the thickness of the last scattering surface (c.f. below), 
the dark energy and dark matter content of the Universe (to an unprecedented accuracy) 
(c.f. figure confirming the earlier claims from supernovae la data and provided 
evidence for early reionization {z ~ 20), which, at least from the point of view of large 
scale structure formation, excludes Warm Dark Matter particle theory models. 

In this review we shall first describe briefly the above-mentioned measurements, and 
then use them to constrain certain particle physics supersymmetric models (in partic- 
ular, the minimally supersymmetric model, constrained by its embedding in a minimal 
supergravity model (mSUGRA)). We shall give a critical discussion on the derived con- 
straints, and discuss the capability of observing Supersymmetry in colliders after these 
latest WMAP data. We shall pay particular attention to discussing regions of the param- 
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eter space of the models, which if reahzed in nature, would imply impossibility of observing 
supersymmetry at LHC and/or linear colliders, even in the context of the mSUGRA mod- 
els. In this respect we shall also discuss the importance of the — 2 experiments of 
the muon gyromagnetic ratio, and how certain results, pointing towards a discrepancy 
between the measured — 2 from the value calculated within the standard model, could 
eliminate the above-mentioned precarious regions of parameter space. 

An important comment we would like to make in this review will concern the dark 
energy component (73 % ) of the Universe. The WMAP measured equation of state for the 
Universe p = wp, with p the pressure and p the energy density, implies —1 < u; < —0.78 
(assuming the lower bound for theoretical reasons, otherwise the upper limit may be 
larger I5j. For comparison we note that w = —1 characterizes a perfect fluid Universe with 
non-zero, positive, cosmological constant. As we shall remark, supergravity quintessence 
models do have this feature w — > —1, and it may well be that by exploiting further 
the data on this dark energy component of the Universe one may arrive at the physically 
correct supergravity model which could constrain the supersymmetric particle physics 
models. 

The outline of this review is as follows: In section two we shall discuss briefly the 
WMAP experiment and its measurements, with emphasis on cosmological parameters 
relevant to particle physics models. In particular, we shall describe the measurements of 
the equation of state of quintessence models, in an attempt to shed light on the nature of 
the dark energy. In section three we shall discuss the basics of inflationary models, which 
seem to be favoured by WMAP, but we shall be critical in our discussion, stressing the 
current inadequacies of the observations. In section four we shall embark on the main point 
of our discussion, namely the connection of the WMAP measurements and Supersymmetry. 
Specifically, we shall discuss various models of SUSY Dark Matter, and argue that only 
Cold Dark Matter (CDM) is favoured by the data, given that Hot Dark Matter (neutrino) is 
excluded in view of the strict upper bounds of neutrino masses imposed by the WMAP, and 
Warm Dark Matter models with a low mass gravitino (less than 10 KeV) are also excluded 
on account of evidence for early (z ~ 20) reionization of the Universe. By concentrating, 
for the purposes of this review, on CDM, made exclusively of neutralinos (viewed as 
the Lightest supersymmetric particle (LSP)) [0], we shall describe, in section five, the 
WMAP-derived constraints on the mSUGRA model, and the prospects for detection of 
supersymmetry within the context of this model in future colliders, as well as direct dark 
matter searches. Also we shall discuss briefly how such constraints, when combined with 
proton decay current lower limits, affect grand unified supersymmetric theories, specifically 
flipped SU(5) models. Finally, in section six we shall review critically the above constraints 
in light of the fact that mSUGRA and all the existing particle physics models do not 
take proper account of the observed dark energy component. In this respect we shall 
discuss supergravity/superstring/brane inspired models with relaxing to zero "vacuum 
energies and argue how such models may be used in the future to constrain particle 
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physics supersymmetric models. The issue is of course still unresolved, as it involves 
the yet unsettled theoretical task of determining the precise mechanism of low-energy 
supersymmetry breaking. Conclusions and outlook, outlining possible future directions, 
will be presented in section seven. 

2 WMAP & (ASTRO)PARTICLE PHYSICS 

2.1 What is WMAP 

The Wilkinson Microwave Anisotropy Probe (WMAP) [3j is a satellite experiment (figure 
dedicated to measure temperature fluctuations of the cosmic microwave background 
(CMB), with unprecedented accuracy, as compared to previous COBE measurements, 
which reaches billionth of a Kelvin degree. 




Figure 2: A graphical picture of the WMAP satellite (courtesy of 
http://map.gsfc.nasa.gov/ |. 



The cosmic microwave background is the afterglow radiation left over from the hot Big 
Bang. Its temperature is extremely uniform all over the sky. However, tiny temperature 
variations or fluctuations (at the part per million level) can offer great insight into the 
origin, evolution, and content of the Universe. The light that is reaching us has been 
stretched out as the Universe has stretched, so light that was once beyond gamma rays is 
now reaching us in the form of microwaves (longer wavelength) (figure We can only 
observe light that comes from the time of last scattering, in a completely analogous way 
as in a cloudy sky we can only see light coming from the cloud surface that is facing 
the earth ground (c.f. figure 4). We remind the reader that the existence of the last 
scattering surface is due to the fact that the Universe expands rapidly and, as a result, 
there is significant dilution in the density of matter particles, including radiation. Hence, 
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the scattering rate of matter and also of light is diminished, and eventually stops (freezing 
out), implying a surface of last scattering. 




Figure 3: The expansion of the Universe causes stretching of the wavelengths of light, so 
light emitted at a gamma ray region of wavelengths reaches us today as microwave radiation 
(courtesy of http://map.gsfc.nasa.gov/). 

In its first year of running, WMAP already provides a much more detailed picture of the 
temperature fluctuations than its COBE predecessor (figure ISJ, which can be analyzed to 
provide best fit models for cosmology, leading to severe constraints on the energy content 
of various model Universes, useful for particle physics, and in particular super symmetric 
searches. 

Theoretically [7j, the temperature fluctuations in the CMB radiation are attributed 
to: (i) our velocity w.r.t cosmic rest frame, (ii) gravitational potential fluctuations on 
the last scattering surface (Sachs- Wolf effect), (iii) Radiation field fluctuations on the 
last scattering surface, (iv) velocity of the last scattering surface, and (v) damping of 
anisotropies if Universe reionizes after decoupling. 

The fluctuations due to (i) imply a dipole anisotropy, while (ii)are the dominant effects 
for large angular scales 9 >> 1°. The fluctuations (iii) - (v) on the other hand, are 
dominant for 6 « 1°. The value oi 6 = 1° is determined by the physics of decoupling. 

The CMB anisotropies in the Sky are expanded as: 

^j, +00 +£ 

1=2 m=-e 

where Yiynid,(p) are spherical harmonics; isotropy implies < aim >= 0. Let us denote the 
variance of aim by =< |a£mP >• A non zero Ci would imply temperature fluctuations. 
More accurately, within a Gaussian model with width a, the above sum over i in the 
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PRESENT 
13.7 Bitlion Years 
after the Big Bang 

The cosmic mkrowave background Radiation's 
"surface of iast scatter" is analogous to the 
iiglvt coming througl^ tlie clouds to our 
eye on a cloudy day. 




We can only see 
the surface of the 
cloud where light 
was last scattered 



Figure 4: We can only see the light emitted from the surface where it was last scattered. 
This happens when we look either at a cloudy sky, or at the cosmic microwave background 
radiation ( courtesy of http://map.gsfc.nasa. gov^ . 



expression for 6T/T is effectively cutoff at I ^ a: 

6T{xi) 6T{x2) 1 , ,s I ,2 -(i±l)2^2 

Ca =< ^ ^ >= ^ 2^{2£ + 1) < \ai,r^r > Pe{xi,X2)e ^ 2 ) ^ 

1=2 

This Gaussian model of fluctuations is in very good agreement with the recent WMAP 
data (see figure inj. The perfect fit of the first few peaks to the data allows a precise 
determination of the total density of the Universe, which as we shall discuss next implies 
its spatial flatness. 

2.2 WMAP Measurement of Cosmological Parameters 

The measurements of the WMAP "F^ on the cosmological parameters of interest to us here 
can be summarized in the figuresdlHl see Refs. jHI, with fij = pij pc. Pi is the matter-energy 
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Figure 5: CMB Temperature fluctuations provided by WMAP and comparison with COBE 
(left), Importance of Detail (right) (courtesy of http://map.gsfc.nasa.gov/ ) . 

Angular scale (deg) 
90 2 0.5 0.2 

7000 E I I I I 3 




Multipole moment I 

Figure 6: Red points (larger errors) are previous measurements. Black points (smaller 
errors) are WMAP measurements (Hinshaw, G. et al. {arXiv:astro-ph/0302217]) . 
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density component, pc = 3H'^ /SttGn (with Gj\[ the Newton constant) the critical density 
for the closure of the Universe, determined within the Friedman- Robertson- Walker frame- 
work, and h the rescaled Hubble parameter, defined through H = lOO/i Km s~^ Mpc~^. 
The WMAP satellite can measure directly Q/i^, which is Hubble parameter independent 
by construction, and can also separately measure h. For our purposes in this article we 
concentrate on the measurements of the density of matter and dark energy (see figure IS)), 
as well as on the equation of state of the dark energy component (c.f. figure IHI)- 
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Figure 7: Likelihood Functions for various cosmological parameters measured by 
WMAP J^. 

We observe from figure IHl that the WMAP results constrain severely the equation 
of state p = wp {p =pressure), pointing towards w < —0.78, if one fits the data with 
the assumption —l<w (we note for comparison that in the scenarios advocating the 
existence of a cosmological constant one has w = —1). Many quintessence models can 
easily satisfy the criterion —1 < w < —0.78, especially the super symmetric ones, which 
we shall comment upon later in the article. Thus, at present, the available data are 
not sufficient to distinguish the cosmological constant model from quintessence (or more 
generally from relaxation models of the vacuum energy) . 

The results of the WMAP analysis (alone), including directly measurable and derived 
quantities, are summarized in the tables appearing in figures ITOIl II 

One therefore obtains the chart for the energy and matter content of our Universe 
depicted in figure^ This chart is in perfect agreement with direct evidence on acceleration 
of the Universe (and hence cosmological constant) from Supernovae la Data 1^ (c.f. figure 
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Figure 8: Energy density and matter density content of the Universe as measured by the 
WMAP satellite alone, as well as in combination with other experiments 



12). 

It should be stressed that the interpretation of the supernovae data is based on a best 
fit Friedmann-Robertson- Walker (FRW) Universe 1 : 



O.SS^Af - 0.6S7a ^ -0.2 ± 0.1 , for 9.m < 1-5 



(1) 



with i^M,A corresponding to the matter and cosmological matter densities. Assuming a 
flat model (k=0), r^totai = 1, as supported by the CMB data, the SNIa data alone imply: 



cFlat 



0.281°;°^ {la stat.)l',f. 



[identified syst. 



(2) 



^ , where a is the cosmic scale factor, receives 



The deceleration parameter defined as q ; 
the following form if we omit the contribution of photons which is very small, 

q = ^Qm -^a^ -0.57 < , (l^A ~ 0.7) . 



(3) 



Hence and ((21) provide evidence for a current era acceleration of the Universe. At this 
stage it should be stressed that the recent observation of two more supernovae at 2; > 1 [21 
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Figure 9: Equation of State of the Dark Energy component of the Universe as measured 
by the WMAP satellite alone, as well as in combination with other experiments 

supports the geometrical interpretation on the existence of a dark energy component of 
the Universe, and argues rather against the role of nuclear (evolution) or intergalactic dust 
effects. 

The recent data of WMAP satellite lead to a new determination of ritotai = 1.02 it 0.02, 
where f^totai = Ptotax/Pc-, due to high precision measurements of secondary (two more) 
acoustic peaks as compared with previous CMB measurements (c.f. figure IHI). Essentially 
the value of Vl is determined by the position of the first acoustic peak in a Gaussian 
model, whose reliability increases significantly by the discovery of secondary peaks and 
their excellent fit with the Gaussian model [S]. 



3 WMAP AND INFLATION 

3.1 Basics of Inflationary Models 

The recent determination of the cosmological parameters by the WMAP team |5j favours, 
by means of best fit procedure, spatially flat inflationary models of the Universe [S]. In 
terms of Robertson- Walker metrics 

ds^ = _de + a^i^t)(-^^^+r^d^A , (4) 
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Table 1 

Power Law ACDM Model Parameters- WMAP Data Only 



Parameter Mean (68% confidence range) Maximum Likelihood 



Baryon Density 




0.024 ± 0.001 


0.023 


Matter Density 




0.14 ±0.02 


0.15 


Hubble Constant 


h 


0.72 ± 0.05 


0.68 


Amplitude 


A 


0.9 ±0.1 


0.80 


Optical Depth 


T 


0.166l°o:S^? 


0.11 


Spectral Index 


Us 


0.99 ± 0.04 


0.97 

1431/1342 



"Fit to WMAP data only 



Figure 10: Cosmological parameters measured by WMAP (only): directly measurable quan- 
tities l^. 

such models are described by an early phase in which the scale factor a{t) is exponen- 
tially expanding with the Robertson- Walker time t: a{t) ~ e^^, where H is the Hub- 
ble parameter. Recent WMAP data confirm its value during inflation to be of order 
H = const = W~^Mpi, where Mpi is the Planck energy scale (~ 10^^ GeV). Spatially 
flat models are described by A; = 0, since the spatial curvature is proportional to this 
parameter, and they correspond to a total energy density f^totai = 1- 

Such a spatially flat Universe is preferred in inflationary scenaria j9j. An important 
feature of conventional inflationary models is the presence of a scalar mode, the inflaton 
4> ^0], which couples to a de Sitter Universe. The field's slow rolling down the hill of 
its potential (c.f. figure IT^ causes the Universe to expand exponentially. Its equation of 
motion in a de Sitter background has a frictional form, with the friction provided by the 
Hubble parameter H (assumed almost constant during inflation): 

^ + 3/70+^ = 0. (5) 
d(p 

where the overdot denotes derivative w.r.t Robertson- Walker time t. 

The quantum fluctuations of a Planck-size (10~^^ meters) Universe, which is assumed 
to be the state of the Universe right after the big bang in inflationary scenaria are enor- 
mously enhanced (c.f. figure IT^ by the exponential expansion and result to cosmic mi- 
crowave background radiation anisotropies. They may also provide the explanation of the 
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Table 1 

Derived Cosmological Parameters 



Parameter 


Mean (68% confidence range) 


Amplitude of Galaxy Fluctuations 


0-8 = 


-- 0.9 ±0.1 


Characteristic Amplitude of Velocity Fluctuations 


o-gfi 


= 0.44 ±0.10 


Baryon Density/Critical Density 


fii = 


= 0.047 ±0.006 


Matter Density/Critical Density 




= 0.29 ±0.07 


Age of the Universe 


to = 


: 13.4 ±0.3 Gyr 


Redshift of Reionization (ionization fraction, = 1) 




= 17±5 


Redshift at Decoupling 


^dec 


= 10881^ 


Age of the Universe at Decoupling 


idee 


= 372 ± 14 kyr 


Thickness of Surface of Last Scatter 


Azaec = 194 ± 2 


Thickness of Surface of Last Scatter 


Atdcc = 115 ± 5 kyr 


Redshift at Matter/Radiation Equality 


^eq ' 


= 3454±|f, 


Sound Horizon at Decoupling 


'"s = 


: 144±4 Mpc 


Angular Diameter Distance to the Decoupling Surface 


dA = 


= 13.7 ±0.5 Gpc 


Acoustic Angular Scale(Z^ = Trd^/rs) 


Ia -- 


= 299 ± 2 


Current Density of Baryons 


ni, = 


= (2.7 ±0.1) X lO"'' cm"^ 


Baryon/Photon Ratio 


V = 


(6.5±°:^) X 10-10 



Figure 11: Cosmological parameters measured by WMAP (only). -derived quantities 0/. 
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0.01 0.10 1.00 

Figure 12: The only Direct evidence for acceleration of the Universe [and hence Dark 
energy) so far is provided by high redshift (z < 1 ) supernovae measurements (SNIa) JIJ/, 
which were prior to WMAP data. 

presently observed large scale structure of the Universe. 

The (thermal) history of the inflationary Universe is quite different from the one 
expected in standard cosmology (c.f. figure E}- During the exponential expansion the 
Universe cools down significantly, and then is reheated after the end of inflation. The 
physics of reheating is not well understood at present, and many scenaria exist in the 
literature on this issue, some of which involve fundamental strings. In the present article 
we would not like to discuss such issues, as they lie outside our scope. An important 
characteristic of inflation is the enormous entropy production at the end of inflation. 

3.2 WMAP-Measured Inflationary Parameters 

The measurements of the WMAP satellite which support inflation are associated with the 
so called Running Spectral Index, a quantity defined as follows |llj . 
Let one consider density fluctuations in the Universe, 

6{x) = 

P 

where barred quantities denote average over the Universe volume V. Taking the Fourier 
transform: 6{x) = j^rj J d^kS^^e^^'^ , one can define the power spectrum of the fluctuations 
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Figure 13: Inflation Basics: the inflaton and its rolling down the hill phase which causes 
inflation (exponential expansion); due to this expansion small quantum fluctuations may 
enhance and be responsible for large scale structures of the Universe at present. 



as [7] 

< 6ix + r)5{x) C = -^V I d'k\6^fe~''^-'^ (6) 

Then, in terms of a comoving scale k, one can define the primordial power spectrum of 
scalar density fluctuations (11) : 

P{k) ^ \6^f . 
The running scalar spectral index ns{k) is 

dXnP{k) 



ns{k) 



dink 



One may fix and its slope w.r.t. k at ko = O.OBMpc ^ and it is customary to assume 
that = 0, from which it follows that 

Then the power spectrum P{k) above is written as 

P{ko) = P{ko) i^- 
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10"*^ lO--"* 10"^'' t(sec) 



Figure 14: The thermal history of the inflationary Universe as compared to the standard 
cosmology: the exponential expansion causes the Universe to cool down significantly, be- 
fore is reheated again. The physics of reheating is not well understood at present. The 
inflationary phase results in an enormous entropy production at the end of inflation. 



The quantity ns{k) carries information on the field derivatives of the inflaton potential 
in best fit models: V , V" , V" {V = and can be measured by WMAP data. 

In fact, standard infiation implies a constant n^, i.e. scale invariance of the fiuctuation 
spectrum. The current value of the running of the spectral index by the available set of 
WMAP data at present is very small but non-zero, however there is not enough statistics 
as yet to arrive at definite conclusions on the shape of the inflaton potential. 

Below we give a list of the basic WMAP measurements that favour inflation For 
deflniteness we give flrst the conventional slow-roll inflationary parameters [7| IH] 

/ \l/2 

where Mp = ( g„-Gjv ) ~ ^'^ ^ ^^^^ GeV is the four-dimensional Planck mass, wth Gn 
the Newton constant, V{(t)) is the inflaton (0) potential, and the prime indicates derivative 
with respect to (j). WMAP has measured the following (c.f. figs. llUlllf) 013]: 



Flatness of the Universe {k = 0) 



Density fiuctuations Ajj, which are expressed in terms of ^ as: Ajj 
2.95 X lO^^A, WMAP result: A = 0.77 ± 0.07 



2 _ V 



247r2M|,e 
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• WMAP results alone are too crude \12\ to determine the shape and height of inflaton 
potential. When, and only when, the WMAP data are combined with all other data, 
then they seem to exclude, at 3 — o" level, a y ~ 0^ potential [3]. 

• The ratio of tensor to scalar perturbations (at quadrupole scale) r = ^ = 16e, 
WMAP result : r ~ 0.16 ^ 

• The Spectral Index for (i) scalar density perturbations = 1 — 6e + 2r/ WMAP 
result: ~ 0.96, and (ii) for tensor perturbations = — 2e . 

• The running of Scalar Spectral Index, a quite important WMAP result, = 
§[(n, - 1)2 - V] + 2^, WMAP result: ^^8x . 

3.3 WMAP and the Inflaton 



T ,0-'- 



TR„[GeVl 



Trh [GeV] 



Figure 15: Branching ratios of various Lepton- Flavour- Violation processes in scenaria 
viewing the Sneutrino as the Inflaton 1131. 

Despite the agreement, via best fit scenaria, of the WMAP data with inflationary 
models, unfortunately no information can be obtained at this stage on the nature of the 
inflaton field. In fact, all the WMAP data imply that at an early stage of our Universe 
there was an era of exponential expansion, but cannot shed light on what caused this 
expansion. Although the most popular and successful model of inflation to date involves 
the inflaton however this is not the only possibility. An exponential expansion of 
the Universe does not have to be produced by an extra fundamental scalar field. The 
inflationary phase may be dynamical in nature, without the need of inflaton fields. Such 
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issues can only be understood when the fuh non-equilibrium dynamics of inflationary 
Universe is at hand. 

We mention in passing at this point that within the context of string theory jjl] , the 
dilaton can play the role of the inflaton fleld, and that this mode may arise dynamically 
in some non-critical string scenaria by incorporating properly non-equilibrium dynamics 
via the so-called Liouville mode 151. The Liouville mode is a world-sheet field jl6j . 
which arises from the need for restoration of the conformal invariance of the underlying 
stringy a- model theory, which is broken as a result of a cosmically catastrophic event, 
such as an initial quantum fluctuation or, in a modern context, the collision between 
brane Universes, one of which represents our world jl7l 1181 I19j . In such models, the 
Liouville fleld appears with a time-like signature, and hence may be viewed |2n( 121) as 
an (irreversible) target time variable in the system. The irreversibility is guaranteed by 
world-sheet renormalization-group arguments, given that the zero mode of the Liouville 
field plays the role of a world-sheet scale [20] . 

Due to its specific coupling with the two-dimensional curvature of the world-sheet of 
the stringy a-model, the Liouville mode contributes linear (in target time) terms in the 
dilaton excitation of the string spectrum. Linear-in-target-time dilaton cosmology has 
been discussed in 22 , where it was shown that it implies linearly expanding Universes 
in the Robertson- Walker time ^ To obtain infiation in such non-critical string scenaria 
requires more complicated string backgrounds JSl ^1 ^1 and such models have 
been discussed in the literature, also in connection with a current-era acceleration of the 
Universe. We shall come back to such issues in section 6. 

Quite recently, other interesting particle physics scenaria have been suggested jl3j . in 
which the infiaton may be the supersymmetric partner of the neutrino (sneutrino). The 
model is argued to be consistent with the WMAP data, corresponding to a quadratic 
inflaton potential y ~ c/)^. In such scenaria Leptogenesis leads to a prediction for the 
Lepton Flavour Violation (LFV) of the model which appears consistent with the WMAP 
constraint on the baryon density at a 3-a level: 7.8 x 10~^^ < < 1.0 x 10"^*^. It must 
be noted however that the model is not consistent with the WMAP non-zero value of the 
running of the spectral index dng/dlnk. However at present the data are too crude for an 
unambiguous determination of such quantities, and so the authors of [T^ used only the 
la level result of that particular measurement to flt their model. 

One can calculate in this model LFV branching ratios BR{t fij) ^ 10^^, BR{^ — > 
67) ^ 10^^^, -Bi?(r 67) < 10^'^, as well as electron and muon dipole moments, which 
have been argued [lH| to be consistent with present experimental limits. The most sensitive 

^In stringy cosmology, the Robertson- Walker frame is defined as tlie one in which the time in the metric 
has unit coefficient, but also the Einstein scalar curvature term in the low-energy effective field-theory action 
has no exponential dilaton terms in front. This requires some redefinition, not only of the time variable, 
but also of the target space metric background by appropriate exponential dilaton factors !22?. Under such 
redefinitions, an originally linear in target time dilaton, becomes logarithmic in Robertson-Walker time. 
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experiments to measure such BRs are those involving neutrino oscihations. Note that such 
BR cannot be measured in the ATLAS/LHC experiment, whose sensitivity to such LFV 
processes is much poorer {BR{t —>■ /U7)atlas ~ 10~^). 

4 WMAP AND DARK MATTER 

4.1 Hot and Warm Dark Matter Excluded 
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Figure 16: Hot Dark Matter models are excluded in view of the WMAP determination of 
(low) upper bounds of neutrino masses j^. 



The WMAP results on the cosmological parameters discussed previously disfavor 
strongly Hot Dark Matter (neutrinos), as a result of the new determination of the up- 
per bound on neutrino masses. 

The contribution of neutrinos to the energy density of the Universe depends upon the 
sum of the mass of the light neutrino species E] : 



94.0 eV 



(9) 



where the sum includes neutrino species that are light enough to decouple while still 
relativistic. 

The combined results from WMAP and other experiments fS^ on the cumulative like- 
lihood of data as a function of the energy density in neutrinos is given in figure [THl Based 
on this analysis one may conclude that fij^/i^ < 0.0067 (at 95% confidence limit). Adding 
the Lyman a data, the limit weakens slightly 5 : 



Q^h^ < 0.0076 



(10) 



or equivalently (from @): "^irn^- < 0.69 eV, where, we repeat again, the sum includes 
light species of neutrinos. This may then imply an average upper limit on electron neu- 
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trino mass < rrii, >e< 0.23 eV. The upper bound on the relevant densities (fTU)) strongly 
disfavours Hot Dark Matter scenarios. 
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00 02 0.4 0* o. 



Figure 17: Numerical simulations for structure formation in Cold Dark Matter ( CDM) and 
Warm Dark Matter (WDM) models with mass mx = 10 KeV \2^ . Upper two figures: 
The projected gas distribution in the CDM (top) and the WDM (middle) simulations at 
z = 20. Lower bottom figures: the distribution of dark halos with mass greater than IO^Mq 
for the CDM (left) and for the WDM (right) model. 

Another important result of WMAP is the evidence for early re-ionization of the Uni- 
verse at redshifts z ~ 20. If one assumes that structure formation is responsible for 
re-ionization, then such early reionization periods are compatible only for high values of 
the masses mx of Warm Dark Matter . Specifically, one can exclude models with mx < 10 
KeV US] based on numerical simulations of structure formation for such models. Such 
simulations (c.f. figure I17|) imply that dominant structure formation responsible for re- 
ionization, for Warm Dark Matter candidates with mx < 10 KeV, occurs at much smaller 
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z than those observed by WMAP. In view of this, one can therefore exclude popular parti- 
cle physics models employing light gravitinos {mx — 0.5 KeV) as the Warm Dark Matter 
candidate. 

It should be noted at this stage that such structure formation arguments can only place 
a lower bound on the mass of the Warm Dark Matter candidate. The reader should bear 
in mind that Warm Dark Matter with masses mx > 100 KeV becomes indistinguishable 
from Cold Dark Matter, as far as structure formation is concerned. 

4.2 Cold Dark Matter in Supersymmetric Models: Neutralino 

After the exclusion of Hot and Warm Dark Matter, the only type of Dark matter that 
remains consistent with the recent WMAP results |S] is the Cold Dark Matter , which in 
general may consist of axions, superheavy particles (with masses ~ 10"^^^^ GeV) j24| I25j 
and stable supersymmetric partners. 

Indeed, one of the major and rather unexpected predictions of Supersymmetry (SUSY), 
broken at low energies MsusY ~ 0(1 TeV), while ii-parity is conserved, is the existence 
of a stable, neutral particle, the lightest neutralino (x?)) referred to as the lightest su- 
persymmetric particle (LSP) Such particle is an ideal candidate for the Cold Dark 
Matter in the Universe [S] . Such a prediction fits well with the fact that SUSY is not only 
indispensable in constructing consistent string theories, but it also seems unavoidable at 
low energies (~ 1 TeV) if the gauge hierarchy problem is to be resolved. Such a resolution 
provides a measure of the SUSY breaking scale MsusY ~ 0(1 TeV). There is indirect 
evidence for such a low-energy supersymmetry breaking scale, from the unification of the 
gauge couplings [22] and from the apparent lightness of the Higgs boson as determined 
from precise electroweak measurements, mainly at LEP |27| . 

This type of Cold Dark Matter will be our focus from now on, in association with the 
recent results from WMAP on relic densities |H1 [3] . The WMAP results for the baryon and 
matter densities (including dark matter) are given in the table of figure ITTIl After combining 
the WMAP data with other existing data, one obtains the following final values: 

Qmh'^ = 0.135i[{;|][J^ (matter) , 

flbh^ = 0.0224 ± 0.0009 (baryons) . (11) 

Hence, assuming that CDM is given by the difference of these two, one has 

^CDAih^ = 0.1126t[]:0i^} , (2 -a level) . (12) 

As mentioned already, in supersymmetric (SUSY) theories the favorite candidate for 
CDM is the lightest of the Neutralinos Xi (SUSY CDM), which is stable, as being the 
Lightest SUSY particle (LSP) From the above-described WMAP results, then, assuming 

^There are cases where the stau or the sneutrino can be the hghtest supersymmetric particles. These 
cases are not favoured )28) and hence are not considered. 
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^CDM — '^^^ infer stringent limits for the neutralino x relic density: 

0.094 < n^h'^ < 0.129 , {2 -a level) (13) 

It is important to notice that in this inequality only the upper limit is rigorous. The lower 
Limit is optional, given that there might (and probably do) exist other contributions to 
the overall (dark) matter density. 

The constraints in SUSY models we shall review in this work come precisely from the 
use of this constraint, together with other results mainly from the anomalous magnetic 
moment of the muon, 5^ — 2, the decay b — > 57, as well as LEP-induced Higgs mass lower 
limits. Specifically, we shall constrain the minimal SUSY Standard Model (MSSM) with 
universal values for the masses of scalars mo the gauginos mi/2 trilinear soft couplings 
^0 at the GUT scale. Furthemore Electroweak Symmetry Breaking (EWSB) is assumed. 
This is the so-called Constrained MSSM (CMSSM). For the purposes of this review we shall 
concentrate on plots {'mQ,mi/2) for the minimal supergravity model (mSUGRA) |29| LSOj. 
in which such CMSSM are embedded. 

It is imperative to notice that all the constraints we shall discuss in this review are 
highly model dependent and the results cannot be extrapolated to non-minimal models. 
But most of the analysis can be extrapolated to such models, with possibly different results. 
In what follows, therefore, for definiteness we shall concentrate exclusively to mSUGRA 
models. We shall not discuss non-minimal models in this article. 



5 WMAP AND SUSY CDM 

5.1 CMSSM/mSUGRA: Basic Features 

Before embarking into a detailed analysis of the constraints of the minimal supersymmetric 
standard model embedded in a minimal supergravity model (CMSSM) [^, we consider 
it useful to outline the basic features of these models, which will be used in this review. 

The embedding of SUSY models into the minimal supergravity (mSUGRA) model 
implies that there are five independent parameters: Three of them, the scalar and gaugino 
masses mo, mi/2 as well as the trilinear soft coupling Aq =, at the unification scale, set the 
size of the Supersymmetry breaking scale. In addition one can consider as input parameter 
tan/3 = ^^^^ , the ratio of the v.e.v's of the Higgses H2 and Hi giving masses to up and 
down quarks respectively. The sign ( signature) of the Higgsino mixing parameter n is 
also an input but not its size which is determined from the Higgs potential minimization 
condition 

U H = P h — 4 — ■ (14) 

^ 2 tan/32-1 tan/32-1 ^ ^ 
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In (|14|) Mz is actually the running Z - boson mass and niHi 2 the soft masses of the Higgses 
H2 and Hi respectively. Si^2 on the r.h.s of this equation are contributions arising from 
loop corrections of the effective potential PTI 1^ . 

The parameter space of mSUGRA can be effectively described in terms of two branches 
as can be seen by exploring the minimization condition p4j) . Following the discussion 
in [23] we can demonstrate this analytically ignoring, for simplicity, the coupling of the b 
quark. In this case p4jl can be approximated by |S21 '■ 



Ciml + C^m% + C'^Al + A/x^p = + ^^1 (15) 



where 



and 



= my, + \aS^ , C'2 = C,-\^ (16) 



C3 ' ' ^4^3 



1 32^0-1 , e 

(1 ^*)' ^2-^ 



A/iioop denote the loop correction to the parameter fi, t = tan/3, Si^2 were defined before, 
the functions e, /, g, k are as defined in [25) and = 1 — {mt/mj)'^ where nif ^ 200 sin/3 
GeV. 

• Ellipsoidal Branch (EB) of Radiative Symmetry Breaking : 

For small to moderate values of tan/3 ^ 7, the loop corrections are typically small, 
and the renormalization group analysis shows that both C2 > and C3 > 0. In this 
case one finds Ci > independently of any scale choice Q for the radiative elec- 
troweak symmetry breaking. Then, one finds that the radiative symmetry breaking 
constraint demands that the allowed set of soft parameters mo and m'1/2 he, for a 
given value of ;U, on the surface of an Ellipsoid. This places upper bounds on the 
sparticle masses for a given value of $ = /x^/A/^ + 1/4. This is the ellipsoidal branch 
(EB) of the Electroweak Radiative Symmetry Breaking (EWRSB). 

Although this is always true for values tan/3 ^ 7, it does not necessarily hold for 
larger values tan/3 ^ 7. In the latter case mo and vn! may lie on the surface 
of a Hyperboloid. This situation, which peresents phenomenological interest, is 
described in more detail below [221 ■ 

• Hyperbolic Branch (HB) of Radiative Symmetry Breaking: 

For large values of tan/3 ^ 7, the loop corrections to [i are significant and [i also 
exhibits a significant variation with varying the scale Q. If we choose the running 
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scale at Qq where the loop corrections are minimized, or even vanish, the loop 
corrections to the r.h.s of H14|) can be ommited and all quantities appearing in eq. 
((TH) should be considered at this scale. Usually the scale Qo is the geometric average 
of the stop masses which give the largest contributions to 5]i^2- In some regions of the 
parameter space it may occur that the sign of Ci is flipped, i.e. sign(Ci((5o)) = ~1) 
so that the minimization condition receives the form 

:.±1. (18) 



In CHI) a = '^^°^^|c^~^'^°' , 13 = '^"^"^^ic,^'^'^"' ^ and the ± sign is determined by 
the condition sign (($ + i)M| - C2AI) = ±1 . It is apparent from ^ that unlike 
the previous case (ttiq, ^'1/2) '^^^ surface of a hyperboloid and hence the 

name Hyperbolic Branch (HB). 

What is interesting in the HB case it the fact that mo and or mi/2 can become very 
large while much smaller values for fi. A subset of HB is the so called high zone. In this 
case EWSB can occur in regions where mo and mi/2 can be in the several TeV range, with 
much smaller values for the parameter fi which however is much larger than Mz- This has 
important consequences for phenomenology as we shall see. In this zone the lightest of the 
neutralinos, xi^ is almost a Higgsino having mass /i. This is called inversion phenomenon 
since the LSP is a Higgsino rather a Bino. 

Except the high zone where the inversion phenomenon takes place the HB includes the 
so called Focus Point (FP) region |351 136j . which is defined as a region in which some 
renormalization group (RG) trajectories intersect (FP region would be only a point, were 
it not for threshold effects which smear it out). We stress that the FP is not a fixed point 
of the RG. The FP region is a subset of the HB limited to relatively low values of mi/2 
and values of /i close to the electroweak scale, Mz, while mo can be a few TeV but not 
as large as in the high zone due to the constraints imposed by the EWSB condition. The 
LSP neutralino in this region is a mixture of Bino and Higgsino and the Higgsino impurity 
allowes for rapid s-channel LSP annihilations, resulting to low neutralino relic densities at 
experimentally acceptable levels. 

From the previous discussion concerning the HB region it is obvious that the crucial 
property characterizing this region is the inversion of the sign of Ci. For large tan/3 ( 
tan/? ^ 7 ), the sign of Ci is mainly driven by that of the quantity 3Dq — 1 as is seen 
from eq. (|17j) . In order for Ci to be negative 3Dq — 1 should be positive which with the 
approximations adopted so far entails to 

. ^ /s" mf rut 
smp > \ - 



2 200 GeV 163.3 GeV 
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Even from this crude approximation it is apparent the central role the value of the running 
top mass mt it plays for the location of the Hyperbolic Branch. If it happens rrit > 
163.3 GeV, at the scale Qq, the HB is never reached. Thus the preference towards smaller 
values of the top quark mass, close to the lowest experimental limit, for the location of 
the HB is manifest from this discussion. Therefore the HB region is sensitive to variations 
of the parameters of the theory which may induce a shift of the theoretically determined 
mt to higher values by 1-2 %. It then becomes evident that two loop RGE effects to 
Yukawa couplings, radiative corrections to the top mass and supersymmetric corrections 
to the bottom Yukawa coupling, which are important for large values of tan/3, as well as 
low energy threshold effects to the strong coupling constant, which greatly affect the RGE 
running of the parameters of the theory and hence low energy predictions, must be duly 
taken into account for a correct and unambiguous phenomenological description of this 
region. These subtleties may drastically affect the conclusions reached in various analyses 
and this is the reason this phenomenologically interesting region is sometimes considered 
as being fine tuned. 

5.2 Focusing The Cosmologically Relevant Regions 

Having described the two principal regions, EB and HB, with different phenomenolological 
characteristics each, we now pass to briefly outline the regions which posses special interest 
as far as CDM is concerned. 

In the EB region the EWSB condition drives the parameter n to be of the same order 
of magnitude with mo, 'rni/2- the bulk of this region the lightest of the neutralinos, 
Xi, is mostly a bino except a region with small mi^2 which Xi has a sizeable Higgsino 
component. This region is however experimentally ruled out by 6 ^ 57 data and light 
chargino searches. Excluding this region we focus our attention to the following subregions 
which possess special features as far as Cosmology and particle physics is concerned. The 
first is the stau coannihilation region in which Xi is almost degenerate in mass, but heavier, 
with the light r-slepton. When this happens the calculation of the neutralino CDM 
requires the coannihilations Xi ^ and f — f [HZl OHl OH] • In this region the mass of Xi is 
allowed to be extended up to 500 GeV assuming that the WMAP data are respected. The 
second region is characterized by large tan /5 ^ 45, and concerns the region where 2 m^o ~ 
niA where is the pseudoscalar Higgs mass. In this region s-channel pseudoscalar 
annihilation toa6 — 5orar — f pair becomes important in opening up extended regions 
compatible with WMAP data. Interestingly enough this large tan (3 region yields high 
values, ~ 10~*'^ pb for the elastic Xi ' nucleon cross section close to the sensitivity limits 
of on going and future direct dark matter experiments. We postpone a more detailed 
discussion of this issue later on in this review. 

Moving to the HB, in general the scalar mg and gaugino mi/2 masses can get very 
large, which implies that the squarks and sleptons can be very heavy, with masses that 
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lie in the several TeV range. The feature of the large mo mass is also shared by the focus 
point region of mSUGRA models [351 136j . Consider the high zone of the HB, in which 
mi/2 ^ ^ ^ Mz- In this scenario one finds [SUES] that the two lightest neutralino states, 
Xi, X2 the light chargino state Xi essentially degenerate with mass ~ /i. Theis 
masses range from several hundreds GeV to above TeV while their mass differences are 
much smaller in the range 1-10 GeV. To first order they are degenerate all having a mass 
/i but this degeneracy is lifted by small corrections 0{^^) entailing to a mass spectrum 

m^o < m + < m^o 

Al Xl A2 

with the lightest chargino states and X2 the next to lightest neutralino. As men- 
tioned earlier in this region the inversion phenomenon takes place given that the lightest 
neutralino switches from being mostly a Bino to being a higgsino. 

The inversion phenomenon has dramatic effects on the nature of the particle spec- 
trum and SUSY phenomenology in this HB. Indeed, in mSUGRA one naturally has 
coannihilation with the sleptons when the neutralino mass extends to masses beyond 
150-200 GeV with processes of the type x^r ^ tl,tZ,th, t&, and tj.^^ 

£"'F','y'y,^Z, ZZ,W^W~ ,hh where I is essentially a f. Remarkably the relic density con- 
straints can be satisfied on the hyperbolic branch also by coannihilation. However, on the 
HB the coannihilation is of an entirely different nature as compared with the stau coan- 
nihilations discussed previously: instead of a neutralino-stau co-annihilation, and stau - 
stau in the HB one has co-annihilation processes involving the second lightest neutralino 
and chargino states [101, Xi ~ xf > followed by Xi ~ X2^Xi ~ Xi ^xf ~ X2 ■ Some of the 
dominant processes that contribute to the above coannihilation processes are |40| 

Xixt^X^xt ^ Uid^eiiy,, AW+, ZW+, W+h 

xtxT^XlX2^U^^^,d^di,W+W- (19) 

Since the mass difference between the states xt ^^^d Xi is the smallest the Xixt coanni- 
hilation dominates. In such cases, the masses niQ Tni/2 may be pushed beyond 10 TeV, 
so that squarks and sleptons can get masses up to several TeV, i.e. beyond detectability 
limits of immediate future accelerators such as LHC. 

Except the high zone the other interesting region belonging to the HB is the Focus 
Point (FP). As already mentioned in the previous chapter unlike the high zone this is 
characterized by mo in the few TeV range, low values of mi 1/2 << mo and rather small 
values of close to Mz- The LSP neutralino in this case is a mixture of Bino and Higgsino 
and its Higgsino impurity is adequate to give rize to rapid s-channel LSP annihilations 
so that the neutralino relic density is kept low at experimentally acceptable values. Since 
fi is small the lightest chargino may be lighter than 500 GeV and the FP region may be 
accessible to future TeV scale colliders. Also due to the relative smallnes of mi/2 in this 
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region gluino pair production may occur at a high rate making the FP region accessible 
at LHC energies. 

As we remarked before although the HB may be viewed as fine tuned, nevertheless 
recent studies |41| . based on a analysis, have indicated that the WMAP data, when 
combined with data on 6 — > 57 and — seem to favour the Focus Point HB region and 
the large tan /3 neutralino resonance annihilation of mSUGRA. 

5.3 SUSY CDM in the pre- WMAP era and the importance of the 
muon's anomaly 

Before embarking onto a discussion on the constraints on SUSY CDM implied by the recent 
WMAP data on relic densities, we consider it as instructive to review phenomenological 
works before the publication of these data. 

The issue of the Dark Matter has attracted the interest of many physicists and in 
the literature there are numerous theoretical works dealing with this issue [421 143L I44L 1451 
Eni mi EHl 113 EDI El HO]. Departures from the minimal mSUGRA models and studies 
including non-universalities IHJ, CP - violating effects IHH as well as efforts related to 
Yukawa unification j381 153j have been pursued. Phenomenological studies including the 
effects of variations of the WIMP velocities on the detection rate have also been considered 
j54j and regions where coannihilation phenomena are important have been analyzed [551 

EOlEZlEHlEniiMlEZI 

The importance of the smallness of the Dark Matter (DM) relic density in constraining 
supersymmetric predictions had been pointed out in the past jSHl EH EOl paying 
special attention at the to the large tan/3 regime. In this region the neutralino (x?) 
pair annihilation through s-channel pseudo-scalar Higgs boson (A) exchange, leads to 
an enhanced annihilation cross sections reducing significantly the relic density [021 • The 
importance of this mechanism, in conjunction with the cosmological data which favour 
small values of the DM relic density, has been stressed in |i58i i59j . 

In previous analyses regarding DM direct searches 121] , we had stressed that the contri- 
bution of the CP-even Higgs bosons exchange to the LSP-nucleon scattering cross sections 
increases with tan p. Therefore in the large tan /3 region one obtains the highest possible 
rates for the direct DM searches and the smallest LSP relic densities. Similar results are 
presented in Ref. IHSl- Early studies on small Q^^h'^ ^ 0.18, taking into account constraints 
from the anomalous magnetic moment of the muon, have been carried out in j,61 , based 
on the mSUGRA model. The results from such analyses are depicted in figure ITHl 

The pre- WMAP situation involving the efforts of other groups, concerning prospects 
for detecting SUSY is summarized in figure TH^ where we also give some benchmark sce- 
naria p^ . 

Undoubtedly one of the most significant experimental results of the last years is the 
measurement of the anomalous magnetic moment of the muon |66| . Deviation of its mea- 
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Figure 18: mSUGRA/CMSSM constraints before WMAP E^/. The Mass of the top 
quark is taken 175 GeV. Dark Green shaded areas correspond to neutralino relic densities 
in the range: 0.08 < ^x^'^ ^ 0.18, whilst light green areas correspond to 0.18 < VL^K^ < 
0.30. The Solid Red Line indicates constraints implied by the anomalous magnetic moment 
of the muon a^^^ = (306ibl06) x 10^"*^^ . The Dashed red line indicates the boundary of the 
region for which the lower bound is moved to 2a limit. The Dashed Blue lines correspond 
to constraints implied by the Higgs mass 113.5GeV < rriHiggs < 117.0GeV. Finally Cyan 
shaded region on the right are excluded due to b ^ sj constraint. 



sured value from the Standard Model (SM) predictions is evidence for new physics with 
Supersymmetry being the prominent candidate to play that role. Adopting Supersymme- 
try as the most natural extension of the SM, such deviations may be explained and impose 
at the same time severe constraints on the predictions of the available SUSY models by 
putting upper bounds on sparticle masses. Therefore knowledge of the value g^ — 2 
is of paramount importance for Supersymmetry and in particular for the fate of models 
including heavy sparticles in their mass spectrum, as for instance those belonging to the 
Hyperbolic Branch. 

Unfortunately the situation concerning the anomalous magnetric moment is not clear as 
some theoretical uncertainties remain unsettled as yet. At present, there are two theoreti- 
cal estimates for the difference of the experimentally measured value of = {g^ — 2)/2 
from the theoretically calculated one within the SM [67j . 

• Estimate (I) a^f^ - af^ = 1.7(14.2) x lO^^^ [0.4(15.5) x 10^^°] 

• Estimate (II) a^f^ - af^ = 24.1(14.0) x 10-^° [22.8(15.3) x 10-^°] 

In (I) the T-decay data are used in conjuction with Current Algebra while in (II) the 
e~e^ — > Hadrons data are used in order to extract the photon vacuum polarization which 
enters into the calculation of g^ — 2. Within square brackets are the most updated values 
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Figure 19: Left figure: SUSY detection prospects before WMAP data l()4h explorable at 
LHC using missing energy and jet signature. Right Figure: Post-LEP Benchmark SUSY 
scenaria Benchmark points have f^^^^ larger than the bounds implies by the WMAP 
data. 



of Ref. [nZI ^- Estimate (I) is considered less reliable since it carries additional systematic 
uncertainties and for this reason in many studies only the Estimate (II) is adopted. Esti- 
mate (II) includes the contributions of additional scalar mesons not taken into account in 
previous calculations. 

In order to get an idea of how important the data on the muon anomaly might be we 
quote Ref. ^33! where both estimates have been used. If Estimate (II) is used at a l.Sfi 
range much of the HB and all of the inversion region can be eliminated. In that case the 
usually explored region of SUSY in the EB is the only one that survives, which, as we 
shall discuss below, can be severely constrained by means of the recent WMAP data. 

On the other hand, Estimate (I), essentially implies no difference from the SM value, 
and hence, if adopted, leaves the HB, and hence its high zone (inversion region), intact. 
In such a case, SUSY may not be detectable at colliders, at least in the context of the 
mSUGRA model, but may be detectable in some direct dark matter searches, to which 
we shall turn to later in the article. 

For the above reasons, it is therefore imperative to determine unambiguously the muon 
anomalous magnetic moment 5^ — 2 by reducing the errors in the leading order hadronic 

■^Due to the rapid updates concerning Qi^ ~ 2 the values of a^f — used in previous works quoted in 
this article may differ from those appearing above. 
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contribution, experimentally, and improving the theoretical computations within the stan- 
dard model. In view of its importance for SUSY searches, it should also be necessary to 
have further experiments in the future, that could provide independent checks of the mea- 
sured muon magnetic moment by the E821 experiment |66j . 

5.4 WMAP and SUSY constraints 
5.4.1 Constraints in the EB 

After the first year of running of the WMAP satellite |3j , the situation concerning mSUGRA 
model constraints in the EB of radiative symmetry breaking parameter space, has changed 
dramatically |68[ I69j . The constraints on the neutralino relic density (|12[ I13|) , implied by 
the best fit Universe model of the WMAP data combined with other experiments [SI E] , 
and with other constraints such as the post-LEP Higgs mass limits, 6 ^ 57 and — 2, 
have restricted severely the available parameter space of the EB. 
Some important remarks are in order at this stage: 

In Ref. inni in calculating the Xi relic abundance, the authors solve the Boltzmann 
equation numerically using the machinery outlined in Ref. jHH]- In this calculation the 
coannihilation effects, in regions where approaches in mass the LSP, which is a high 
purity Bino, are properly taken into account. 

The details of the procedure in calculating the spectrum of the CMSSM can be found 
elsewhere |7U1 1611 IBU] , Here we shall only briefly refer to some subtleties which turn out 
to be essential for a correct determination of the Higgsino, /i, and Higgs, m|, mixing 
parameters, which greatly affect the Neutralino, Chargino and Higgs mass spectra, espe- 
cially in the region of large tan/3. This region is of particular phenomenological interest as 
we discussed in previous sections, /i and mg are determined by minimizing the one-loop 
corrected effective potential (HJ. For large tan/3 the derivatives of the effective potential 
with respect the Higgs fields, which enter into the minimization conditions, are plagued by 
terms which are large and hence potentially dangerous, making the perturbative treatment 
untrustworthy. In order to minimize the large tan/3 corrections we had better calculate 
the effective potential using as reference scale the average stop scale ~ ^/^ii • 
At this scale these terms are small and hence perturbatively trusted. 

The proceedure in solving the renormalization group equations for masses and cou- 
plings follows essentially that of Ref. [72 . The low energy threshold corrections to the 
strong coupling constant are duly taken into account and the unification scale Mqut is 
defined at the point where the gauge couplings ai,2 meet. No unification of the gauge cou- 
pling constant is enforced with 01^2 at Mqut as this usually results to too high values 
for the strong coupling constant at Mz- Therefore to comply with the experimentally 
value for as{Mz) we take it as input and thus the value of as at the unification scale is 
an output rather than an input. 
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A significant correction, which drastically affects the numerical proceedure arises from 
the gluino-sbottom and chargino-stop corrections to the bottom quark Yukawa coupling 
hb [731 1741 rm I75j . The proper resummation of these corrections is important for a correct 
determination of hh [7^1 EZi and this has been taken into account. 

Seeking a precise determination of the Higgs boson mass the dominant two-loop correc- 
tions to this have been included [75]. Concerning the calculation of the 6^57 branching 
ratio, the important contributions beyond the leading order, especially for large tan P, 
have been taken into account (75] . 

There have been two independent groups working on this update of the CMSSM in 
light of the WMAP data, with similar results j681 169] and below we summarize the results 
of their analysis in figures 1201211 for some typical values of the parameters tan/3 and 
signature of /u. 

As mentioned in previous sections the region designated as large tan/3 has particular 
phenomenological interest. In this region the dominant neutralino annihilation mechanism 
is through s - channel annihilation of the pseudoscalar Higgs boson (A), XiXi ~^ i'b or tt. 
The reason is that by increasing tan /? the mass tjia decreases, while the neutralino mass 
remains almost constant, if the other parameters are kept fixed. Thus the pseudoscalar 
Higgs mass ttia is expected eventually to enter into the pole regime in which = 
2m^o. At the same time the pseudoscalar Higgs couplings with rr, bb are large as being 
proportional to tan/3. Thus although the x? is a high purity bino, which does not couple 
to A, its small Higgsino component, which does couple, can give sizeable effects due to 
the largeness of the Abb and Arf couplings. Both effects result to enhanced annihilation 
cross sections dominating the process and hence to small neutralino relic densities |69] over 
a broad region of the parameter space characterized by large values tan/3 values. From 
this discussion it becomes obvious that for a correct determination of the relic density, 
in the large tan/3 region, an unambiguous and reliable determination of the A-mass, niA, 
is required. In the CMSSM, niA is not a free parameter but is determined once the 
other parameters are given. niA depends sensitively on the Higgs mixing parameter, m^, 
which is determined from the minimization conditions of the effective potential, and for 
its calculation all 1 - loop corrections must be taken into account. It is found (Hill that 
the handy calculation of the pseudoscalar mass through the second derivatives of the 
effective potential is scale independent and approaches the pole mass to better than 2% 
provided that one includes the corrections of all particles, especially those of Charginos 
and Neutralinos, in addition to the leading corrections of the third generation sfermions. 

The analysis of is concentrated on large values of tan/3 for reasons mentioned 
above. It is found (see the right panel of the lower part of figure [2T]) that for such values, 
a region opens up within which the relic density is cosmologically allowed. This is due 
to the pair annihilation of the neutralinos through the pseudo-scalar Higgs exchange in 
the s-channel. In this region the ratio m^/2m approaches unity and the pseudo-scalar 
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exchange dominates as explained above. It is for this reason that we give special emphasis 
to this particular mechanism which opens up for large tan f3 and delineates cosmologically 
allowed domains of small relic densities and large elastic neutralino - nucleon cross sections. 
In this case the lower bound put by the (g/j, — 2) data cuts the cosmologically allowed region 
which would otherwise allow for very large values of mo,mi/2- 
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Figure 20: mSUGRA/CMSSM constraints after WMAP 1681. Dark Blue shaded region 
favoured by WMAP ( 0.094 < il^/i^ < 0.129 ). Turquoise shaded regions have 0.1 < 
^x^^ — Brick red shaded regions are excluded because LSP is charged. Dark green 
regions are excluded by b ^ 57. The Pink shaded region includes 2 — a effects of g^ — 2. 
Finally, the dash-dotted line represents the LEP constraint on e mass. 
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Figure 21: Cosmologically allowed regions of the relic density for in the {mi/2,mQ) plane for 
large tanj3 values '69L The mass of the top is taken 175 GeV. In the dark green shaded area 
0.094 < il^o hi < 0.129. In the light green shaded area 0.129 < il^o < 0.180 . The solid red 
lines mark the region within which the super symmetric contribution to the anomalous magnetic 
moment of the muon is a^^^^ — (361 ± 106) x 10^^^. The dashed red line is the boundary of 
the region for which the lower bound is moved to its 2a limit. The dashed-dotted blue lines are the 
boundaries of the region 113.5 GeV < rnuiggs < 117.0 GeV. The cyan shaded region is excluded 
due to b ^ s"f constraint. 
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In figure !^ we see that for the tan (3 = 55 case, close to the highest possible value, and 
considering the 2 a lower bound on the muon's anomalous magnetic moment a^'^^^ > 
149 X 10~^^ and values of fl^o in the range 0.1126^Q'Qj^gj^ , allowed points are within 
a narrow stripe. The point with the highest value for mo is (in GeV) at (1710,1711/2) = 
(850,550) and that with the highest M1/2 at(mo,mi/2) = (750,600 ). The latter 
marks the lower end of the line segment of the boundary 149 < 10^^^ a^^"^^ which 
amputates the cosmologically allowed stripe. 

It should be noted that within la of the E821 data only a few points survive which lie 
in a small region centered at (mo, rni/2) = (725, 300) . The bounds on mo, mi/2 displayed 
in (the lower part of) figure TTH refer to the j4o = case. Allowing for Aq ^ values, the 
upper bounds put on mo,Mi/2 increase a little and so do the corresponding bounds on 
sparticle masses. For the LSP, the lightest of the charginos, stops, staus and Higgses the 
upper bounds on their masses are displayed in Tabled for various values of the parameter 
tan /? , if the new WMAP determination [3 of the Cold Dark Matter (O and the 
2a bound 149 < 10"^^ a^^^^ < 573 of E821 is respected. We have also taken into 
account the limits arising from Higgs boson searches as well as from 6^57 experimental 
constraints. In extracting these values we used |^ a random sample of 40, 000 points in 
the region \Aq\ < 1 TeV, tan/3 < 55, mo,Mi/2 < 1-5 TeV and /i > . The lightest of 
the charginos has a mass whose upper bound is ~ 550 GeV , and this is smaller than the 
upper bounds put on the masses of the lightest of the other charged sparticles, namely the 
stau and stop, as is evident from Tabled Hence the prospects of discovering CMSSM at a 
e^e~ collider with center of mass energy ^/s = 800 GeV, are not guaranteed. Thus a center 
of mass energy of at least ^/s 1.1 TeV is required to discover SUSY through chargino 
pair production. Note that in the allowed regions the next to the lightest neutralino, x' , 
has a mass very close to the lightest of the charginos and hence the process e^e~ xx' ^ 
with x' subsequently decaying to x + or x + 2 jets, is kinematically allowed for such 
large tan /3, provided the energy is increased to at least ^/s = 860 GeV. It should be noted 
however that this channel proceeds via the t-channel exchange of a selectron and it is 
suppressed due to the heaviness of the exchanged sfermion. Therefore only if the center of 
mass energy is increased to y/s = 1.1 TeV supersymmetry can be discovered in a e~^e~ 
collider provided it is based on the Constrained scenario |69| . 

An important conclusion, therefore, which can be inferred by inspecting the figures 
1201 1211 is that the constraints implied by a possible discrepancy of — 2 from the SM 
value ( a^'^^'^ ^ 15.0 x 10~^° ), when combined with the WMAP restrictions on CDM 
(neutralino) relic densities H13|) . imply severe restrictions on the available parameter space 
of the EB and lower significantly the upper bounds on the allowed neutralino masses m^o. 
This is clearly seen in figure [22 In this the width of the strips is much smaller than the 
spacing between them. For tan/3 < 40 one obtains jHE] n^xi ^ GeV (see figure l23]) . 
which is significantly lower than the bounds obtained in pre- WMAP cosmological models. 



34 



1500 




100 1000 2000 2500 

mi,2 (GeV) 



Figure 22: Strips correspond to tanf3 = 5,10,15,20,25,30,35,40,45,50,55. Light Green 
Regions of the graph {mo,mi/2) are those compatible with the WMAP result on the neu- 
tralino relic density 0.094 < Q^h^ < 0.129. Darker Blue shaded parts of strips are com- 
patible with — 2 results (at a 2 — a level) \()^ . 
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Table 1: Upper bounds, in GeV, on the masses of the lightest of the neutralinos, charginos, staus, 
stops and Higgs bosons for various values of tan/? if the new WMAP value 0^ for flcDAih^ o.'^d 
the 2a E821 bound, 149 x 10"" < a;^^'^^ < 573 x 10~", is imposed \U^ . 

To summarize, therefore, in the EB the recent WMAP data [3j can easily be adapted 
through by most of the benchmark points (in bulk and co- annihilation regions) of figure 
El by reducing mp. In the post- WMAP benchmark scenaria 81 the only exception is 
the point H (tip of co- annihilation tail): WMAP data will lower thereby facilitating 

detection prospects of SUSY at LHC or future linear e+e~ colliders a center of mass energy 
y/s = 1.1. TeV |69(l68j. provided the HB region, and in particular its high zone (inversion), 
is not realized in nature. 
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Figure 23: The Neutralino Mass Range in the EB of radiative symmetry breaking parameter 
space, implied by the WMAP data and other constraints (b — > 37, f;^ — 2 etc), in the context of 
the mSUGRA model \6i^ . Left: < 0, Right: n > 0. Red solid lines: Upper limits without — 2. 
Dashed Blue lines: Upper limits with g^ — 2. Dotted lines: Green: Masses of , Blue masses of 
cl, at the tip of co- annihilation tails. 



5.4.2 Constraints in the HB 

Despite the above-mentioned good prospects of discovering minimal SUSY models at fu- 
ture colliders, if the EB is realized, however, things may not be that simple in Nature. A 
recent study |1J of mSUGRA in light of the recent WMAP data has indicated that the 
HB/focus point region of the model's parameter space seems to be favoured along with the 
neutralino resonance annihilation region for /x > and large tan/3 values. The favoured 
focus point region corresponds to moderate to large values of the Higgs parameter fj,'^, and 
large scalar masses mo in the several TeV range. The situation is summarized briefly in 
figure 1211 

The situation in case the HB is included in the analysis is depicted in figure HJl '3^ , 
where we plot the mo — rni/2 graphs, as well as graphs of mo, mi/2 vs the neutralino LSP 
mass. The neutralino density is that of the WMAP data. We stress again that, in case the 
high zone (inversion) region of the HB is realized, then the detection prospects of SUSY 
at LHC are diminished significantly, in view of the fact that in such regions slepton masses 
may lie in the several TeV range (see figure 0^. 

In view of the above results, an updated reach of LHC in view of the recent WMAP 
and other constraints discussed above (see figure 0^)1 has been performed in |52], showing 
that a major part of the HB, but certainly not its high zone, can be accessible at LHC. 
The conclusion from this study jS2] is that for an integrated luminosity of 100 fb~^ values 
of mi/2 ~ 1400 GeV can be probed for small scalar masses mg, corresponding to gluino 
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Figure 24: WMAP data seem to favour (-^ < 4/3J (green) the HB/focus point region 
(moderate to large values of fj,, large niQ scalar masses) for almost all tan(3 (Left), as well 
as s - channel Higgs resonance annihilation (Right) for /_i > and large tan(3 



masses rUg ~ 3 TeV. For large mo, in the hyperbolic branch/focus point region, mi/2 ~ 700 
GeV can be probed, corresponding to rUg ~ 1800 GeV. It is also concluded that the LHC 
(CERN) can probe the entire stau co-annihilation region and most of the heavy Higgs 
annihilation funnel allowed by WMAP data, except for some range of ttt-q, mx/2 iii the case 
tan/3 ^ 50. 

A similar updated reach study in light of the new WMAP data has also been done for 
the Tevatron ,83!, extending previous analyses to large rriQ masses up to 3.5 TeV, in order 
to probe the HB/focus region favoured by the WMAP data [3Tj. Such a study indicated 
that for a 5(7 (3(t) signal with 10 (25) fb~^ of integrated luminosity, the Tevatron reach 
in the trilepton channel extends up to mi/2 ~ 190 (270) GeV independent of tan/3, which 
corresponds to a reach in terms of gluino mass of mg ~ 575(750) GeV. 

5.5 WMAP and Direct Dark Matter Searches 

We turn now to study the impact of the new WMAP data on Dark Matter, the constraints 
from {g^ — 2) and b ^ s'y, and the Higgs mass bounds on direct Dark Matter searches. A 
quantity of great interest is the spin-independent ( scalar ) neutralino-proton cross section 
on which experimental limits exist from the current dark matter experiments (T^0p(5'I) < 
10~^ picobarns. 

In figure |27| we plot the scalar Xi-nucleon cross section as function of the LSP mass, 
m^o inn]- On the top of it the shaded region (in cyan colour) is excluded by the CDMS 
experiment ISH. The DAMA sensitivity region is plotted in yellow 1851. Pluses (+) (in blue 
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Figure 25: mo — graph, and tuq and mi/2 vs. graphs, including the HB of 

mSUGRA ^J. Such regions are favoured by the WMAP data piy . 
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(in red 



colour) represent the rest of the points of our random sample which consists of 40,000 
points in the range |^o| < 1 TeV,mo,mi/2 < 1-5 TeV, tan/3 < 55 and /i > belonging 
to the EB region. Here the Higgs boson mass, rrih > 113.5 GeV and the b ^ sj bounds 
have been properly taken into account [Hni- From this figure it is seen that the points 
which are compatible both with the (g^ — 2) E821 and the cosmological data can yield 
cross sections slightly above 10"^ pb when m^o is about 120 GeV. The maximum value of 
m^o is around 200 GeV but in this case the scalar cross section drops by almost an order 
of magnitude 10~^ pb. Accepting the 2 a (g^ — 2) bound the maximum value of the scalar 
cross section is again ~ 10~^ pb, for m^o ~ 120 GeV, but the m^o bound is increased 

Al Al 

to about 280 GeV at the expense of having cross sections slightly smaller than 10~^ pb. 
Considering the /i > case, it is very important that using all available data Q^, one can 
put a lower bound « 10~^ pbarns on the scalar cross section which is very encouraging for 
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mSiigra with Um^ - VO, A^j - 0, ^ > 




Figure 26: Left: The updated Reach in {1^10,1711/2) parameter plane of mSUGRA assuming 
100 fb^^ integrated luminosity. Red (magenta) regions are excluded by theoretical (ex- 
perimental) constraints 182L Right: Contours (in view of the uncertainties) of several 
low energy observables : CDM relic density (green), contour of nih = 114.1 GeV (red), 
contours o/a^lO"'^'^ (blue) and contours of b ^ BF (xl^"^) (magenta). 

future DM direct detection experiments ISS]. Such a lower bound cannot be imposed when 
fi < 0, since the scalar cross section can become very small due to accidental cancellations 
between the sfermion and Higgs exchange processes. However, this case is not favoured 
by iOfi — 2) and b ^ s^ data. 

Portions of the high zone (inversion) region of the HB may be detected by such direct 
Dark Matter Search Experiments, which are thus complementary to collider physics. As 
in the EB case such scattered plots can be produced also for the HB jSH], which seems 
to be favoured by the WMAP data [IJ, and which has regions that lie beyond the reach 
of LHC. As stressed in [,VA\ it is important to distinguish the spin-dependent from the 
spin-independent (scalar) cross sections in this branch. 

In figure I2E1 we give a plot of the scalar cross section for tan/3 = 10 and /i > 0. We do 
observe that the lower rapidly falling curve that terminates at = 300 GeV is the branch 
on which staus co-annihilation occurs. The upper curve arises from the low zone of the HB, 
while the patch to the right is the one corresponding to the inversion region of the HB |33j . 
This patch indicates that the scalar cross sections are quite significant even though one is 
in the inversion region in the range 10"^'' — 10~^ pbarns. Thus, despite the fact that the 
direct detection of SUSY in the inversion region is difficult, the neutralino-proton scalar 
cross section are still substantial. The sensitivity of the future dark matter detectors [HHl 
will reach 10~^ pbarns, being capable of probing a significant portion of the parameter 
space of fig. I28L In this latter respect we also mention novel indirect detectors of SUSY 
dark matter exploiting possible detector of cosmic antimatter (antideuteron) , as a 
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Figure 27: Neutralino-nucleon scalar cross section vs. in the EB 16!^ . 

result of neutralino interactions (see figure I29|) , which could detect neutralinos in a mass 
range from 80-350 GeV, much higher than other future dark matter experiments j86j . 

Finally, in figure EOl we plot the spin dependent parts of the neutralino-proton cross 
section, for the same range of parameters as in figure 1281 A comparison between these 
two figures shows that the spin dependent cross section is much larger, by three or four 
orders of magnitude, than the spin independent one. One can repeat the analysis to find 
that for tan/3 < 50 the neutralino mass range consistent with the WMAP constraints on 
the branch corresponding to neutralino-stau coannihilation is m-^o < 500 GeV |681 169j , 
and m^o < 1200 GeV for the high zone of the HB, where the relic density constraints are 
satisfied due to coannihilation with the next to lightest neutralino and light chargino , 
as discussed previously. We repeat again that these constraints remain intact under the 
imposition of the Estimate (I) of the g'^ — 2, while the constraint arising from the inversion 
region of the HB is removed under the imposition of the Estimate (II) of the — We re- 
iterate, therefore, once more the importance of having unambiguous measurements of this 
quantity in the immediate future, which will be indispensable in guiding super symmetric 
searches in future colliders. 

5.6 WMAP, Proton Decay and SUSY Constraints in Grand Unified 
Supersymmetric Models 

In view of the above-described constraints on MSSM embedded in mSUGRA, it is also 
natural to ask whether additional constraints are implied by the lower limit on the proton 
life time which is a prediction of Grand Unified Theories (GUTs). So far we did not 
consider grand unification in the above analysis. 

There are a number of important issues associated with such an extension. One is 
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Figure 28: Spin independent part of the neutralino-proton cross section vs. m^o including 
the HB W^ . 

the exact value of sin^ 6w, which acquires important corrections from threshold effects at 
the electroweak scale, associated with the spectrum of MSSM particles [HSllEn]) and at 
the grand unification scale, associated with the spectrum of GUT supermultiplets |H51inn] . 
Precision measurements indicate a small deviation of sin^ 6\y even from the value predicted 
in a minimal supersymmetric SU{5) GUT, assuming the range of as{Mz) now indicated 
by experiment |91j . 

The second issue is the lifetime of the proton. Minimal supersymmetric SU (5) avoids 
the catastrophically rapid p — > e~^TT^ decay that scuppered non-supersymmetric SU{5). 
However, supersymmetric SU (5) predicts p vK^ decay through d = 5 operators P^PIH] 
at a rate that may be too fast to satisfy the presently available lower limit on the 
lifetime for this decay [031001 • The latter requires the SU{5) colour-triplet Higgs particles 
to weigh > 7.6 x lO^^GeV, whereas conventional SU{5) unification for as{Mz) = 0.1185 it 
0.002, sin'^9w = 0.23117±0.00016 and aem{Mz) = 1/(127.943±0.027) [SEl would impose 
the upper limit of 3.6 x 10^^ GeV at the 90% confidence level |94| . This problem becomes 
particularly acute if the sparticle spectrum is relatively light, as would be indicated if the 
present experimental and theoretical central values of {g^ — 2) (SHI remain unchanged as 
the errors are reduced. 

The simplest way to avoid these potential pitfalls is to flip SU{5) [03 EH]- As is well 
known, flipped SU{5) offers the possibility of decoupling somewhat the scales at which the 
Standard Model SU{3), SU{2) and U{1) factors are unifled. This would allow the strength 
of the U{1) gauge to become smaller than in minimal supersymmetric SU{5), for the same 
value of as{Mz) [HOj- Moreover, in addition to having a longer p — > ej ix^if' lifetime than 
non-supersymmetric ^[/(S), flipped 5f7(5) also suppresses the d = 5 operators that are 
dangerous in minimal supersymmetric SU{5), by virtue of its economical missing-partner 
mechanism j97j . 
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Figure 29: GAPS Detector of Indirect SUSY Dark Matter Searches, through the detection 
of cosmic antideuteron, may cover a neutralino mass range that may be complementary to 
all other direct dark matter and collider searches 187}/ . 

In [M] , the authors re-analyzed the issues of sin^ 6w and proton decay in flipped 
SU{5) IHH) in view of the most recent precise measurements of as{Mz) and sin^ 6w, 
and the latest (at the time of publication) limits on supersymmetric particles (before 
publication of WMAP data). The analysis of such issues was done in the framework 
of CMSSM, as above, and consisted of both a general analysis in the {mi/2,mo) plane 
and also more detailed specific analyses of benchmark CMSSM parameter choices that 
respected all the available experimental constraints at the time |65| . It was found that 
the p — > e/fi'^TT^ decay lifetime exceeds the present experimental lower limit [HS]) with 
a significant likelihood that it may be accessible to the next round of experiments jlUUj . 
In the context of the present article we shall not present the analysis in great detail, in 
particular we shall not discuss the ambiguities and characteristic ratios of proton decay 
modes in flipped SU{5). We refer the reader to the existing literature for this [221 ■ Instead 
we shall give a comprehensive description of their results in figures 1311 1321 which are self 
explanatory. 

In figure |^ we give the pre- WMAP results of [HHl ) while in figure OH we give the 
updated post WMAP figures (taking into account the updated benchmark points |81j). 
for comparison. From figs.l3H and|S21it is evident that the 'bulk' regions of the parameter 
space preferred by astrophysics and cosmology, which occur at relatively small values of 
{mi/2,rno), generally correspond to r(p —>■ e"'"7r°) ~ (1 — 2) x 10^^ y. However, these 
'bulk' regions are generally disfavoured by the experimental lower limit on the Higgs mass, 
ruh, and/or by 6 — > 57 decay. Larger values of r(p — > e'^Tr^) are found in the 'tail' regions 
of the cosmological parameter space, which occur at large mx/2 where x~ ^ coannihilation 
may be important, and at larger mi/2 and mo where resonant direct-channel annihilation 
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Figure 30: Spin dependent part of the neutralino-proton cross section vs. m^o including 
the HB U^ . 

via the heavier Higgs bosons A, H may be important. 

Another issue studied in [S9j were the possible imphcations of the GUT threshold 
effect (5heavy |HH1 EOl- A general expression for this in flipped SU{5) is given in [55] . 
where we refer the interested reader. In fig. |2S1 we display [H^ the possible numerical 
effects of ^heavy on t{jp e/fi'^TT^) in the various benchmark scenarios, assuming the 
plausible ranges —0.0016 < (5heavy < 0.0005 B9:. The boundary between the different 
shadings for each strip corresponds to the case where (5heavy = 0- The left (red) parts 
of the strips show how much r(p e'^vr'^) could be reduced by a judicious choice of 
Shenvy, and the right (blue) parts of the strips show how much r(p — > e~^7r^) could be 
increased. The inner bars correspond to the uncertainty in sin^ 6w- On the optimistic 
side, we see that some models could yield t{p e'^ir^) < 10^^ y, and all models might 
have t{p e~^7r^) < 5 x 10^^ y. However, on the pessimistic side, in no model can we 
exclude the possibility that r(p e'^'vr'^) > 10^^ y. 

We notice at this stage that a new generation of massive water- Cerenkov detectors 
weighing up to 10^ tonnes is being proposed |100j . that may be sensitive to r(p — > e^Tr^) < 
10^^ y. According to the calculations of such an experiment has a chance of detecting 
proton decay in flipped SU{5), though nothing can of course be guaranteed. We also 
mention |1031I104] that flipped SU{5) makes predictions for ratios of decay rates involving 
strange particles, neutrinos and charged leptons that differ characteristically from those 
of conventional SU{5). Comparing the rates for e^, /i"*" and neutrino modes would give 
novel insights into GUTs as well as mixing patterns. 

We conclude from the above analyses therefore, that flipped SU{5) evades two of 
the pitfalls of conventional supersymmetric SU{5): (i) it offers the possibility of lowering 
the prediction for as{Mz) for any given value of sin^ 6w and choice of sparticle spec- 
trum, and (ii), as far as proton decay is concerned, flipped SU{5) suppresses p — > 
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Figure 31: Pre-WMAP constraints on supersymmetric flipped SU(5) model ^99]. The solid (blue) 
lines are contours of t{p e/zi+Tr*^) in the {mi^2Ti^a) plane for the CMSSM with (a) tan/3 — 
10,^ > 0, (b) tan/3 = 10,^ < 0, (c) tan/3 = 35, ,u < and (d) tan/3 = 50,^ > 0. The (blue) 
crosses indicate the CMSSM benchmark points with the corresponding value of tan /3 and sign of 
H f6'5y . Following 11 (Jl^ . the dark (red) shaded regions are excluded because the LSP is charged, the 
light (turquoise) shaded regions have 0.1 < ft^^h^ < 0.3, intermediate (green) shaded regions at low 
™i/2 o-re excluded byh ^ sj, shaded (pink) regions at large (mx/2j"^o) are consistent with — 2 at 
the 2-a level, and electroweak symmetry breaking is not possible in the hatched regions. The near- 
vertical dashed (black) lines correspond to the LEP lower limit on chargino mass, = 103.5 GeV, 
the dot-dashed (red) lines to the Higgs mass, ruh = 114 GeV as calculated using the FeynHiggs 
code \10S}/ . and the dotted (blue) lines at small 't^o) to mg = 100 GeV. 
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Figure 32: As in fiaure \!^l\ but taking into account the new cosmological constraints from WMAP 
on relic densities (dark blue regions). The contours of t{p e/zi+Tr*^) are given here by the 
light-blue solid lines (we are grateful to J. Walker for providing us with these updated figures). 
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Figure 33: For each of the CMSSM benchmark points, this plot shows, by the lighter outer bars, 
the range of t{p e/ii^TT^) attained by varying (5hoavy over the range -0.0016 to + 0.0005 ]ES- 
The central boundary of the narrow inner bars ( red, blue ) corresponds to the effect of ^ught alone, 
with ^hoavy = 0, while the narrow bars themselves represent uncertainty in siv? 9w ■ We see that 
heavy threshold effects could make t{p ^ e//i+7r*') slightly shorter or considerably longer. 

decay naturally via its economical missing-partner mechanism. As in conventional su- 
persymmetric SU{5), the lifetime for p e/n'^n^ decay generally exceeds the present 
experimental lower limit. However, as shown in 99 , the flipped SU{5) mechanism for 
reducing as{Mz) reduces the scale at which colour SU{3) and electroweak SU{2) are uni- 
fied, bringing r(p — > e/ ii^tt^) tantalizingly close to the prospective sensitivity of the next 
round of experiments. Proton decay has historically been an embarrassment for minimal 
SU{5) GUTs, first in their non-super symmetric guise and more recently in their minimal 
super symmetric version. The answer may be to flip SU (5) out of trouble. We remark also 
that flipped SU{5) grand unified supersymmetric theories arise naturally in intersecting 
brane models |lU5j . 

6 Dark Energy: Cosmological Constant or 
Quintessence or...? 

Before ending this review we would like to devote a few pages on an important matter, 
which so far has not been properly taken into account. In our constrained minimal models 
above, we have taken into consideration the WMAP results concerning dark matter density 
constraints, which constitutes only 23 % of the Uninerse energy density content. However, 
as mentioned in section two, a plethora of recent data, including those of WMAP, has 
indicated that 73 % of the Universe energy consists of an unknown substance, termed 
dark energy (c.f. figures [Tl8|) . 

The origin of this component is as yet unknown, despite enormous theoretical effort. 
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In the simplest approach, the dark energy is nothing other but a cosmological constant. 
A negative cosmological constant was introduced by Einstein in the original theory of 
relativity to enforce compatibility of Einstein's equation of General Relativity with a static 
Universe, believed (wrongly) to be the case of the observed Universe at the time. Hubble's 
observation on the expansion of the Universe (1927) removed the necessity for introducing 
such a constant into the theory. Ever since, until five years ago, theoretical physicists 
believed that the vacuum energy density of the Universe was zero, and they were trying 
(alas in vain!) to find a symmetry reason for that. The situation changed drastically 
recently, after experimental evidence for a present-era acceleration of the Universe. A 
FRW Universe with a positive cosmological constant constitutes the best fit model of the 
WMAP and supernovae data From a theoretical viewpoint, however, this poses 
a serious threat in quantization. Indeed, a de Sitter Universe, which would be the final 
state of such a positive cosmological constant FRW Universe, due to the expansion of the 
Universe which leads to a dilution of any matter density, is characterized by the presence of 
an event horizon, making the proper definition of asymptotic states, and hence a scattering 
matrix (S-matrix), problematic. This presents serious challenges for a proper quantization 
of the theory, which are unresolved at present jlU6j . It also poses a serious challenge for 
string theory jl()7| 110^ . given that the latter, at least in its initial perturbative form, has 
been formulated as an S-matrix theory. 

Perhaps a more satisfactory, from this point of view, approach to the problem of a 
physical explanation of the currently observed acceleration of the Universe PP, and the 
deduced dark energy component lEJ is to assume a time varying vacuum energy, A{t) 
relaxing to zero asymptotically in the Robertson- Walker time jlU91 111(7] , t — > oo. 

One may obtain this time dependent "vacuum energy" by means of a (scalar) field, 
which has not relaxed yet to the minimum of its potential. Such a mechanism, known 
as quintessence, has attracted a lot of attention |lllj . and represents an elegant non- 
equilibrium approach to the issue of a positive cosmological constant. The important ad- 
vantage of quintessence (or other relaxation) models is the possibility of defining asymp- 
totic states, and hence an 5-matrix, properly, due to the asymptotic vanishing of the 
vacuum energy, when the system reaches its equilibrium. Another important advantage 
is that there exist attractor solutions in a cosmological setting for some classes of poten- 
tials |112j . which imply that in some cases the scalar field will join the attractor solution 
before the present epoch for a wide range of initial conditions, thereby evading (in prin- 
ciple) fine tuning and tracking problems. There are, however, serious issues that remain 
unresolved, most of them associated with the nature of the tracking quintessence field, 
which should be important today, but not in the past of our Universe. Another important 
issue is the current coincidence in order of magnitude of the dark energy density with the 
matter energy density of the Universe, which quintessence models are called to explain in 
a fine-tuning free way. 

A natural relaxation mechanism, compatible with the current value of A(tpi.oscnt) ^ 
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10^^^'^ Mp inferred from the observations, would be the one in which |1U91 IllUj A(t) ~ 
Mp/t^, for long times Mpt S> 1, where Mp ~ 10^^ GeV is the four-dimensional Planck 
energy scale. For Mptpresent ^ 10^*^, which is the Age of the (observed) Universe, one 
obtains naturally the above limits. 

Such inverse square law relaxation models exist in the literature, in either field or 
string/brane quintessence theory jlU9l IllUj . or even in non- equilibrium approaches to 
string theory, such as non-critical Liouville strings jl5[ I17| I18L I19j . where the vacuum- 
energy relaxation law is obtained from first principles, specifically, by means of world-sheet 
conformal field theory methods. 

In this last respect the quintessence field is actually provided by the Liouville mode 
itself, which is viewed as the time variable. It is the nature of non-critical strings to 
associate the Liouville mode with a linear in time part of the dilaton, and hence non- 
critical string quintessence may actually be viewed as a dilaton quintessence, which also 
characterizes some other approaches within critical string theory jll3j . The advantage of 
dilaton quintessence is that the scalar dilaton field exist in the string spectrum, and there 
is no need to assume extra fundamental scalars, as is the case in quintessence models of 
local field theory. 

There is also another aspect of non-critical strings, which makes them interesting, at 
least from a theoretical viewpoint. Unlike critical strings, which are formulated as a theory 
of target-space scattering amplitudes (S-matrix), Liouville strings do not necessarily have 
a well-defined target-space S-matrix j2Ul llU8j . Upon the identification of target-time with 
the Liouville mode |2f)j , the world-sheet Liouville correlators among vertex operators may 
admit the interpretation of non-factorizable $-matrix elements, acting on density matrices 
rather than target-space state vectors ^, pout =^Pm, with S , with S the ordinary 
S-matrix, which may thus not be well defined. However, we stress that the operator $ is 
well defined in Liouville strings, as being associated with correlation functions of vertex 
operators PUI . 

Such non-critical string models are then capable of describing consistently even quan- 
tum cosmological models with a non-zero cosmological constant, which, as we mentioned 
earlier, are hamprered in a field theory context by the lack of an S-matrix formalism 
due to the existence of cosmic horizons. The lack of a proper ^-matrix may also have 
implications on violations of CPT symmetry in such systems jll41 l2Uj , with interesting 
phenomenological consequences in both particle physics and cosmology |115j . Of course, 
in Liouville-quintessence relaxation string models one may be able of defining an S-matrix 

*The presence of density matrices implies an analogy of Liouville strings with open quantum mechanical 
systems |2U| . The analogy is exact, given that in such models there are modes which are not accessible 
by low-energy scattering experiments, such as topological gravitational modes or back-reaction space- 
time foam effects, or, in colliding brane world scenaria, recoil modes etc.. Intgegrating out such modes 
from a low-energy observer, defines an effective low-energy string theory, which is non critical and out of 
equilibrium, ft is such theories that cosmological models are based upon in this context |15lll7llTSlll9| . 
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asymptotically in time t, given the possibility of defining properly asymptotic states at 
the equilibrium point t — > cxd, where critical string theory is reached in such cases. 

As we have mentioned previously, in section 3, in some string theory models the dilaton 
has also been viewed as the inflaton [21 HH ^1 El IS]- It would be nice to find a 
consistent realistic model in which a single field can explain both phases of the Universe. 
From the point of view of the acceleration of the Universe, implied by the dark energy, 
such an assumption is not unnatural given that currently the Universe appears to enter 
an accelerating de Sitter phase which is qualitatively similar to that of (the beginning of) 
inflation. This may be the result of a particular dilaton potential in a string/brane model. 
However, this issue is still far from being resolved, given that it is associated with non 
perturbative effects in string theory which are not at hand. 

Prom the phenomenological point of view quintessence models are still compatible 
with the current observations, including those of WMAP |116[ IH]. Indeed, as discussed 
in section 2 (c.f. figure 1^1), the WMAP best fit model analysis has provided us with the 
following range of the values of the quintessence field Q equation of state: 

^PQ _ {Qf/2-V{Q) 
PQ iQ)V2 + ViQy 
-l<w< -0.78 (20) 

where V{Q) is the potential of the quintessence, and the lower limit is theoretical, corre- 
sponding to the cosmological constant model (if this theoretical assumption is relaxed the 
upper limit becomes slightly larger, ~ —0.67 jH])- 

Such values appear naturally in many quintessence models jllll 1171 1181 I19j . From 
our point of view in this work we would like to concentrate and review approaches that 
embed such quintessence scenaria in supergravity models jll71lll^lll9j . We note that, in 
this conext, the first phenomenological approach was given in |110j . well before the recent 
astrophysical data Such models could in principle be used, then, to constraint 
minimal super symmetric models, such as MSSM, in the way mSUGRA was used above, 
but this time taking into account the dark energy component of the Universe energy 
density ^. 

The issue of looking at quintessence in supersymmetric theories is motivated clearly by 
the need to reconcile particle physics models with cosmology. The presence of cosmological 
constant or dark vacuum energy immediately brings up the issue of compatibility with the 
hierarchy between supersymmetry breaking scale and the size of the vacuum energy density 
of the Universe. In supersymmetric quintessence models such issues may find a natural 
explanation, given that scalar cosmological fields in cosmology, such as quintessence fields 

^The cosmological constant problem in mSUGRA models is not solved. Even if one fine tunes the tree 
level contributions to be small, however loop corrections result in Planck size contributions to vacuum 
energy. 
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and inflatons, may arise naturally in a super symmetric particle spectrum and it is not 
unreasonable to expect that the current height of the relevant potential may become 
compatible with the phenomenologically correct size of supersymmetry breaking (of at 
most a few TeV). 

Scalar field models of quintessence typically require that the expectation value of the 
scalar field today is of order of the Planck Mass, in order to explain the observed accel- 
eration of the Universe. This immediately implies that one should include supergravity 
corrections to such models (in the case of supersymmetric cosmological models). Brax 
and Martin jll7j were the first to discuss supergravity corrections in the quintessence 
potential in a supergravity model. They considered a (toy) supergravity model with a 
superpotential 

VF(Q) = A°+3q-" (21) 

where A is the (mass) scale of the particle theory ^, and a flat Kahler potential K = QQ*, 
where Q from now on denotes a scalar superfield. 

The resulting scalar potential for such a model |117j has an F-term of the general 
form 

F(Q) = F2-e^' ^Sk^IT^P (22) 
which in the case of a flat Kahler potential, gives jll7j : 

,,c)^.#='^(«!^-(, + 4Q^.i.'Q^) ,23) 

with 13 = 2(a + 1), and k = SirMp"^ the gravitational constant. Notice that, for a > 
which is the case of supersymmetric models, the potential has non-polynomial terms in the 
quintessence field, which when combined with the exponential term turn out, as we shall 
discuss below, to be crucial for the correct phenomenology of the model in the present 
epoch. 

However, as it stands, the scalar potential (|23|) is ill defined for < Q >= Mp due to 
negative contributions. One way out jll7j would be the imposition of the constraint that 
the superpotential vanishes < >= 0, which eliminates the dangerous negative terms in 
the scalar potential. 

The model has several attractive features then. The presence of non-polynomial terms 
in the scalar potential (|23|) implies compatibility of the model with a current era acceler- 
ation of the Universe. Moreover, due to the presence of the exponential term, the value of 
the equation of state, loq, computed from is pushed towards the value —1 in contrast 
to the usual non supersymmetric case for which it is difficult to go beyond ujq ~ —0.7. For 
Qrn ~ 0.3, as suggested by the current data, the model of |117j predicts ujq ^ —0.82, which 

®For instance, in supersymmetric SU(N) theories A is the scale at which the gauge couphng becomes 
non perturbative |118| . 
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Figure 34: Left: The evolution of the energy density of the quintessence field in the model 
of ref. 1119] versus that of the background (z is the redshift). Right: The evolution of the 
equation of state of the quintessence field in that model. In the figures, the dash-dotted 
line represents the background energy density, the continuous lines correspond to the initial 
condition for the quintessence field Qm/Qm ^ 1; while the dashed lines correspond to the 
initial condition for the quintessence field Qin/Qm ^ 1 (where Qm denotes the value of 
the quintessence field at the minimum of the scalar potential). 

is in agreement with the recent WMAP data implying —1 < wq < —0.78 (see figures IH)), 
as discussed in section 2. 

The condition < >= is possible to realise in the presence of other matter fields but 
in general is a tight restriction. Moreover, such a constraint is not compatible with models 
of supersymmetry breaking. Things therefore might have been easier if this restriction is 
relaxed. In ref. |119j a different Kahler potential has been used than in |117j : 

K =[ln{KQ + KQ*)]^/K^ . (24) 
The resulting scalar potential has the form: 

V = M^ (2x2 + (4a - 7)x + 2(a - l)^) - e[(i-^)'-2"(i-^)] (25) 

X 

where x = {—^kQ)'^/^, and = A^^'^°' k'^^'^°'2'^°' . In this case the minimum of the 
resulting scalar potential is always positive for a > 1.25 without the need for the imposition 
of any constraint on < >. 

The cosmology of that model is given in figure 1341 It is important to notice that 
the presence of field Q denominators in the potential (|25j) accompanying the exponential 
terms are crucial for yielding a negative equation of state w < in the current era, thereby 
implying acceleration of the Universe today, as seems to be predicted by observations ^El- 
Indeed, a solution for the potential involving only exponential and polynomial factors of 
the quintessence fields cannot yield a negative w today ^. After reaching its minimum the 

^We note at this stage that the case with only exponential factors of the quintessence field in the 
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scalar field will mimic the cosmological constant with wq = —1. At present, the WMAP 
data point towards Wcai — > — 1, which is certainly a feature that characterizes supergravity 
models (c.f. figure 151)) . but the big question remains as to whether one can distinguish 
(supersymmetric) quintessence with w ^ —1 from a cosmological constant model. This 
question is associated clearly with the topic of the present review. If one constructs an 
appropriate supergravity model of quintessence which can incorporate consistently the 
dark energy component of the Universe with a broken supersymmetry at TeV scale then 
one may use such models to constrain particle physics supersymmetry searches as done 
above for the MSSM. From such studies one may eventually make predictions that could 
distinguish a quintessence model from a naive cosmological constant one. This remains 
the biggest challenge for the years to come. The reason is the following. 

At this stage it is worth mentioning that the authors of |119) have also considered 
another type of Kahler potentials, 

K = -\[ln{KQ + KQ*)] . (26) 

Such logarithmic potentials are characteristic of no-scale supergravity models |121j . which 
are known to be derived at low energies from string models |122j . as well as from M- 
theory [HH] . 

However, the resulting scalar potential, in terms of canonically normalized fields of the 
form Q = (IukQ) / V^K, reads: 

V{Q) = (^A5+/3^i+/3l^!_^^ e-^f"^^ (27) 

with P = 2a + 1. Positivity of the potential is guaranteed for P > \/3. 

The problem with this model, however, is that, as it stands, the potential ()27() has the 
form of the above-mentioned scaling solution of quintessence |12Uj , according to which one 
cannot have a current era acceleration of the Universe, since one cannot obtain a negative 
equation of state. We note, however, that such a conclusion is valid only if there is one 
quintessence field. As we shall discuss later JJj, in more complicated string models, with 
several field configurations, one has the possibility of obtaining a current era acceleration 
and a negative value of the equation of state parameter w, by appropriately arranging 
couplings and field magnitudes. Thus, it may well be that no-scale SUGRA models can 

potential is discussed in |120| . and in this case the attractor is a 'scaling one'. However these results 
apply in field theory where there is only a single quintessence field, and may not characterise the case of 
multifield configurations in more complex situations, such as string theory. Indeed, as we shall discuss 
below, a negative quintessence equation of state today w < can be obtained for exponential dilaton 
potentials, at least in a non-critical string cosmology framework |15lil7 18 19^, where the vacuum energy 
scale is set by the central charge deficit of the non-critical cr- model |lti|. and expresses physically a departure 
from equilibrium, due to cosmically catastrophic events in our world. 



52 



still provide realistic, phenomenologically acceptable, cosmologies to date. This remains 
to be seen. 

A serious problem for all the supersymmetric quintessence models so far is Super- 
symmetry Breaking (SSB) at a phenomenologically acceptable scale M5, such that M| ~< 
F > ^ (IQio GeV)2 for gravity-mediated cases of SSB, or M| ~< F > ^ (10^ GeV)^ 
for gauge-mediated cases (with F denoting an F-term contribution, see (EH)). In the model 
of |119j this requires a value for T^~<i^>K~^~ 7713/2/^"^, where 777,3/2 is the gravitino 
mass, in order to cancel the F-term contribution and give negligible vacuum energy, in 
agreement with the known limits (or recent observations) on the cosmological constant (we 
remind the reader that, in this formalism, the cosmological constant is given by < F2 >, 
c.f. ()22|) '). Although the relaxation of the < W >= constraint allows for such solutions, 
however it is clear that, in such models, the dynamical cosmological constant provided by 
the quintessence field cannot be the dominant source of SUSY breaking. The authors of 
jll9j modified their superpotential by adding the term 7773/2 on the right hand side of 
(|21j) but then the scalar potential acquires changes 5V containing this constant term, 

5V ~ A3+"7773/2K-" + 7772/2^"^ . (28) 



This leads to a serious disruption of the quintessence potential, due to the constant second 
term. 




Figure 35: Various non-critical string scenaria of quintessence Left: A non-critical (Liouville) 
string scenario for quintessence involving colliding five-brane worlds, compactified on magnetized 
tori In this model the dilaton field today is constant (with vev of order of the Planck scale), 

while it was not constant in inflation era. The dashed lines in the evolution of the world-sheet cen- 
tral charge deficit are conjectural, in order to match the criterion of consistency with nucleosyn- 
thesis. The continuous lines have been calculated using world-sheet (logarithmic) conformal field 
theory methods. Right: The evolution of the energy density of the quintessence field in the non- 
critical string model of Diamandis et al. \njj versus that of the background. In that model, there is a 
time dependent dilaton on the observable brane world, whose time dependence today is logarithmic, 
and causes the appearance of a relaxing to zero effective vacuum energy K{t) ^ MpC^'^ , <i> ^ —hit. 
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This is still a serious unresolved problem on field theoretic quintessence models. It 
may well be that string theory models tackle such problems due to unconventional su- 
persymmetry breaking mechanisms one may devise in that framework, as a result of the 
extra dimensions. For instance, it is advocated in jl24j that one may actually live in a 
four-dimensional world with unbroken supersymmetry, and that the mass splittings be- 
tween the superpartners occur as a result of the obstruction of supersymmetry, through 
some higher dimensional effects in string/brane theory. By 'obstruction' we mean here a 
situation in which the ground state of the system is still supersymmetric, characterised 
by a zero cosmological constant, but the excitations of the system exhibit supersymme- 
try breaking mass splittings. If this is the case, then one does not have to worry about 
supersymmetry breaking terms in the effective SUGRA quintessence potential. However, 
a concrete realisation of these ideas, in the context of realistic models, is still lacking. 
Another uncnventional idea about SSB is advocated in ref. |125j . according to which the 
relation between the supersymmetry breaking scale Ms and the cosmological constant 
< -F^ > (c.f. (|22|) ') may be modified, in such a way that the amount of vacuum en- 
ergy, claimed to be observed today, is responsible for the (phenomenologically acceptable) 
supersymmetry breaking. This could be the case if one had ther following relation: 

Ms = K-^<F>^Ky , (29) 

with P = 1/8 instead of the usually assumed 1/4. In this case, one would not have any 
dangerous F-terms in the SUGRA potential of the order k~^M|, which would disrupt the 
quintessence potential. At present, however, at least in the context of conventional field 
or string theory, such a scenario remains only a conjecture. 

We now remark that alternative ways out of these problems could be provided by 
completely unconventional ways of looking at stringy cosmology, namely the approach to 
quintessence advocated in |171 1181 119j , by means of non-critical string theory jl5j . As we 
shall discuss below, in such a framework one may also obtain in some cases a concrete 
realisation of with P = 1/8. 

According to such scenaria, the quintessence field is provided by the Liouville mode, 
which is the linear in time part of a string theory dilaton. The current acceleration of the 
Universe is due to an excitation of our brane world, as a consequence of a catastrophic 
cosmic event. For instance, in the model of |18[ I19j such a catastrophic event is provided 
by the collision of two brane worlds, moving initially with a small relative velocity u <^ c 
(c.f. figure l35|) . This results in a departure from criticality in the resulting world-sheet 
string theory describing excitations on the brane. In more general situations, such as the 
models of Diamandis et al. |17j . such a departure from equilibrium could be due to an 
initial quantum fluctuation of a single brane world. 

The vacuum energy is relaxing to zero as Mp/t^ for long times after the collision, or 
the initial brane fluctuation, and measures essentially the deviation of the excited string 
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system from equilibrium. This scaling is obtained by (logarithmic) conformal field theory 
methods on the world sheet of the string describing excitations of the (observable) brane 
world. It is worth stressing that in the model of jl7j . the effective quintessence (dilaton) 
potential is of an exponential form V{^) ~ Mpe^^ , with the dilaton scaling with (large) 
time as ~ —Int. Naively, one is tempted to think that such potentials could not cause 
a current-era acceleration of the Universe |12Uj . since this potential characterises the so- 
called 'scaling solution' in cosmology |12()j . However, this conclusion is only true if there 
are no other fields in the problem. As demonstrated in jJT], the string models considered 
there, have a plethora of other fields (moduli and flux fields), whose magnitudes can be 
arranged in such a way so as to yield a quintessence equation of state with w ^ —0.75, and 
thus a current-era acceleration of the Universe, as claimed by the observations El- At 
present, however, one could not identify precise dynamical reasons for such arrangements 
in the magnitude and the couplings of the various field configurations, and hence, at 
present, such models may be considered as fine-tuned. 

In the above-described models of non-critical string quintessence the basic source of 
supersymmetry breaking is independent of the value of the vacuum energy today. For in- 
stance, in Diamandis et al. ^7], one uses type-0 strings which are not supersymmetric due 
to appropriately projecting out superpartners in the string spectrum |12f)j . On the other 
hand, in the models of jl8l I19j . involving five-brane worlds, the supersymmetry breaking 
is due to the compactification of the two extra dimensions on internal manifolds with mag- 
netic fields. The supersymmetry breaking in this case is obtained as a consequence jl27j 
of the mass splittings in the spectrum of string excitations due to the different coupling 
of fermionic and bosonic excitations on the brane world to the internal magnetic fields 
(the associated Nielsen-Olesen instabilities of such scenaria are welcome when the model 
is considered in a cosmological setting). In fact in some cases within the models of jSj it 
is even possible to realise concretely the scenario of |125j . given that one can derive the 
relation H29|) with (3 = 1/8, although we stress such a relation is not the only phenomeno- 
logically consistent solution in these models. As remarked in however, the biggest 
challenge for such non-critical string theory models of quintessence is the nucleosynthesis 
era, which seems to require an almost critical string setting (c.f. the left picture in figure 

ESI. 

Finally, we wish simply to mention that, in the modern context of brane theories, there 
are new interesting approaches to the cosmological constant problem, without, however, 
offering any solution, at present. These explore either the holographic principle, by means 
of a conformal-field-theory (CFT)/de-Sitter-space (positive cosmological constant) corre- 
spondence (in the line of the correspondence between CFT and Anti-de-Sitter Space (with 
negative cosmological constant) of Maldacena tl2H|)! or condensation of tachyons in (non- 
supersymmetric) string theory /brane models ^. Although we think of such issues as very 

*We note that tachyonic instabilities may be useful in cosmology, as triggering inflationary expansion 
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interesting and intellectually challenging, however, we shall not touch upon them here, 
due to lack of space. We refer the reader to the relevant literature for more details |129j . 

This completes our discussion on dark energy and its currently unsolved problems. 
From our brief but hopefully comprehensive discussion in this section, the reader must 
have certainly realized that the problem is far from being solved, however it presents many 
interesting possibilities for future work, both theoretical and experimental. Theoretically, 
the biggest challenge is to find a correct model that accounts for the currently observed 
acceleration of the Universe, compatible with other constraints from particle physics, such 
as supersymmetry breaking. Such a model can then be used to constrain super symmetric 
particle physics models, in a similar spirit to mSUGRA discussed here. Experimentally, 
it is desirable to have more direct evidence on the current acceleration of the Universe by 
having more dedicated supernovae (or even other type) experiments. This will provide the 
only direct evidence for such a phenomenon, which when combined with indirect evidence 
from CMB anisotropies, will give a definite result on the important issue of dark energy. 



7 CONCLUSIONS AND OUTLOOK 

In this article we have given a comprehensive review of the latest constraints on CMSSM, 
especially after the first data released by the WMAP satellite experiment. Specifically, we 
have combined the new WMAP cosmological data |SI on Dark Matter with recent high 
energy physics experimental information including measurements of the anomalous mag- 
netic moment of the muon, from E821 Brookhaven experiment, the 6 — > 57 branching ratio 
and the light Higgs boson mass bound from LEP and we studied the imposed constraints 
on the parameter space of the CMSSM. We have assessed the potential of discovering 
SUSY, if it is based on CMSSM, at future colliders and DM direct search experiments. 
The available parameter space seems to be constrained severely by the combination of 
the above-mentioned data. The use of the new WMAP data in conjunction with the 2 a 
{g^ — 2) bound can guarantee that in LHC but also in a e^e~ collider with center of mass 
energy ^/s ~ 1.1 TeV CMSSM can be discovered provided that the high (inversion) region 
of the HE is not realized. 

However, WMAP data seem to favour HB/focus point region (this was based on a 
analysis of mSUGRA model). Updated reach studies have indicated that, for an integrated 
luminosity 100/6" \ LHC can probe mi/2 ~ 1400 GeV, for small mo (gluino mass rrig ~ 3 
TeV). For large mo the HB/focus point region mi/2 ~ 700 GeV can be probed (gluino 
mass rrig ~ 1800 GeV) via conventional SUSY channels. 

However HB/focus point region appears to extent indefinitely to large mi/2,nT-o values, 
far beyond LHC reach. Such regions can be probed by direct dark mater searches. 

of the Universe in the past, but being completely harmless today |17II129| . 
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The effect of these constraints is also significant for the direct DM searches. For the 
/i > case we found that the minimum value of the spin-independent Xi-nucleon cross 
section attained is of the order of 10~^ pb p^j . 

We have also discussed such constraints, combined with those coming from proton 
life time current lower limts, in the context of grand unified supersymmetric models, 
specifically the flipped SU(5) model, which seems to survive such tight constraints. 

For future directions we point out that one should take seriously into account the 
other constraint coming from recent astrophysical data, that of dark energy. We believe 
it is quintessence, i.e a relaxing to zero (non-equilibrium) field. The WMAP data point 
towards an equation of state of quintessence w = p/ p ^ —1 which is close to that of a 
cosmological constant. This feature is shared by Quintessence SUGRA models |117j . and 
should be explored further within Superstring, Brane models of a relaxing to zero vacuum 
energy, in conjunction with the issue of supersymmetry breaking. 

The important theoretical task for the years to come will be to determine the (phys- 
ically) correct SUGRA/string/brane model to constrain SUSY searches exploiting the 
non-zero Dark Energy component of Universe. Experimentally, it is certain that the post- 
WMAP era, which now begins, will imply that Particle Physics and Astrophysics will 
proceed together for the exciting years to come and provide useful and complementary 
physics input to each other. 
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